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A systematic expansion of the induced inclusive gluon radiation associated with jet pro- 
duction in a dense QCD plasma is derived using a reaction operator formalism. Analytic 
expressions for the transverse momentum and light-cone momentum distributions are de- 
rived to all orders in powers of the gluon opacity of the medium, Na g /A — L/\ g . The 
reaction operator approach also leads to a simple algebraic proof of the "color triviality" 
of single inclusive distributions and to a solvable set of recursion relations. The analytic 
solution generalizes previous continuum solutions (BDMPS) for applications to mesoscopic 
QCD plasmas. The solution is furthermore not restricted to uncorrelated geometries and al- 
lows for evolving screening scales as well as the inclusion of finite kinematic constraints. The 
later is particularly important because below LHC energies the kinematic constraints signif- 
icantly decrease the non-abelian energy loss. Our solution for the inclusive distribution also 
generalizes the finite order exclusive (tagged) distribution case studied previously (GLV1). 
The form of the analytic solution is well suited for numerical implementation in Monte Carlo 
event generators to enable more accurate calculations of jet quenching in ultra-relativistic 
nuclear collisions. Numerical results illustrating the constributions of the first three orders 
in opacity are compared to the "self-quenching" hard radiation intensity. A surprising result 
is that the induced gluon radiation intensity is dominated by the (quadratic in L) first order 
opacity contribution for realistic geometries and jet energies in nuclear collisions. 
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I. INTRODUCTION AND SUMMARY 



One of the expected signatures of the quark-gluon plasma created at RHIC and LHC energies of y/s ~ 
200 AGeV and ^fs~ ~ 1500 AGeV is jet quenching [ ||-[ ||. At the SPS energies of y/s ~ 20 AGeV, on 
the other hand, no quenching of moderate p± < 4 GeV hadrons [ Q was observed even in Pb + Pb.. This 
may be due to the break-down of pQCD at such low momentum scales and the difficulty of disentangling 



non-perturbative multiparticle production effects such as the Cronin and intrinsic kx effects [ 25 . Also 
finite kinematic constraints may limit the energy loss at low jet energies. In addition, at SPS energies jet 
propagation in high density matter is limited to at most a few fm/c due to rapid longitudinal expansion. This 



motivated our previous work (GLV1) [ |26|, |27j to study the problem of energy loss in "thin" plasmas. In that 
work we considered the exclusive (tagged) case of energy loss associated with a fixed number of interactions. 
Rough numerical estimates for nominal 5 GeV jets suggested that at those relative small energies and small 
opacities, the radiative energy loss may indeed be much smaller than predicted by Baier, Dokshitzer, Mueller, 
Peigne and Schiff (BDMPS) [ [HJ |n], g Q and Zhakharov [ |T| Q for asymptotic jet energies. 

Several approaches have been advanced to compute non-abelian energy loss. One approach is aimed at 
treating relatively "thick" targets, which while small compared to the jet coherence length are large enough 
so that many collisions occur in the medium. This leads to a continuum formulation of the problem. In [ 
|Io| , |TT| , |l2| , |l6| an effective 2D Schroedinger equation formulation was introduced, and in [ O, M [l9| |o| 



a path integral formulation of the problem was developed. Another approach [ ^, ^, |26[ 27, p8| , which we 
extend here, aims to address the problem of radiative energy loss in "thin" plasmas, a few mean free paths 
thick, by computing directly the radiation pattern from the finite number of Feynman diagrams for the case 
of a few collisions. 

The advantage of the former is that it can make direct contact with the conventional Landau-Pomeranchuk- 
Migdal effect [ [30[ [H| in QED. It has the disadvantage that the continuum solutions obtained for the case of 
many collisions cannot be directly applied to the case of a few collisions. Also the effects of finite kinematical 
constraints are difficult to include in such approaches. The advantage of the second approach is that the few 
collision case can be computed directly from the finite number of amplitudes involved. The disadvantage is 
that the LPM limit is out of reach, and the work involved in summing diagrams increases exponentially with 
the number of collisions. In the present paper we introduce a new approach to bridge the gap between these 
approaches. Our new approach is based on the construction of a suitable reaction operator, i?„, from which 
recursion relations for the inclusive gluon distribution can be derived and solved analytically at arbitrary 
order n. 

For "thin" QCD plasmas that can be formed in nuclear collisions, the gluon radiation intensity can be 
studied systematically through an expansion in powers of opacity defined by the mean number of collisions 
in the medium 

L Na e i f f 2 da e i(z) / z\ dN a ei R G 

n = — = — — = / dz / a q — -5 p [z, r = - w — — - „ log , (1) 
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where N is the number of targets in the medium of transverse area An opacity expansion in terms of 
the path integral formulation was introduced in [ [l9|, Q . In [ p^l we considered the exclusive tagged 
case where all N target partons interacting with the jet. Here (and []2^]) we extend that calculation to the 
inclusive case where fluctuations of the number of collisions are allowed and only the geometry is fixed. 

For a homogeneous rectangular target of thickness, L, the density is p = N/(LA±), and the mean free path 
is A = l/(pa e i). The opacity is then simply L/X = Na e i/A±. For application to high transverse momentum 
jets propagating through cylindrical nuclear reaction geometries, we can interpret L — 1,2 A 1 ' 3 = R s . For 
more realistic 3+ ID Bjorken and transverse expanding Gaussian cylindrical geometry such as 

At 2 \ /2|x|Ar 
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L is replaced by the equivalent rms Gaussian transverse radius Rg = 0.75 A 1 / 3 fm. The rightmost expression 
in (Q) is obtained by averaging the number of collisions of a transverse jet in this expanding Gaussian cylinder 
with the initial jet production coordinate, Xo, averaged over p(xo, To). Here To is the formation time of the 
plasma. At RHIC energies (y 7 ! ~ 200 AGeV) the expected rapidity density of the gluons is dN/dy ~ 1000 
for A — 200 [ . For a typical elastic gluon-gluon cross section er e ; ~ 2 mb and a plasma formation time 
~ 0.5 fm/c, the opacity is moderately small h < 10. This suggests that neither the thin nor thick plasma 
approximations may apply for applications to nuclear collisions. Therefore, an alternate method is needed 
to handle the intermediate (mesoscopic) case. 

An important simplification may be anticipated due the non-abelian analog [ |[ [To[ ^8) of the Landau- 
Pomeranchuk-Migdal (LPM) effect [ |3(], |3f|. While the number of scatterings is moderately large, the 
radiation intensity angular distribution and total energy loss are controlled by the combined effect of the 
number of scatterings n — L/X and the formation probability pf ~ L/lf — L\k\ 2 / (2xE) of the gluon in 
the medium. Here x is the light-cone momentum fraction and |k| <C xE is the transverse momentum 
of the radiated gluon. Therefore, the opacity expansion may converge more rapidly than one would first 
expect. One of the surprising results derived here and summarized in [|2Sfl is that the inclusive induced gluon 
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radiation intensity is dominated by the (already quadratic in L) first order opacity contribution for realistic 
geometries and jet energies in nuclear collisions. 

In GLV1 [ p7| , we showed in contrast that in the exclusive tagged case, i.e. with all recoiled partons 
measured in coincidence with the jet and gluon, a much more complex and nonlinear non-abelian radiation 
pattern emerges. The tagged recoil partons act as additional color dipole antennas that interfere in a 
nontrivial way with the jet's color dipole antenna. The resulting color algebra cannot be reduced to simple 
powers of the color Casimirs. We found that the gluon number distribution for the case of n s = N tagged 
target partons has the form 



1 + ((! + £)"• (K,d/Z ns , (3) 



where f( n *) depends on two dimensionless ratios K = k 2 //i 2 and £ = A/// = \fi 2 /2xE, and depends only 
weakly on n s . The color dependence enters through a rapidly increasing function of R — Ca/Cr, the ratio 
of the Casimirs in the adjoint (d,A = N 2 — 1 dimensional) gluon representation and the Casimir, Cr, of 
the da dimensional jet representation. In GLV1 unitarity was imposed by a wave function renormalization 
factor, Z ns , that was computed perturbatively. For the inclusive case studied here, unitarity is assured by 



the inclusion of certain contact double Born (Virtual) terms as emphasized in BDMS [ 12 

To address energy loss of hard probes produced in nucleus-nucleus collisions we concentrate as in [ [27 , 28 
on the case of hard jet produced inside the plasma at a finite point (io,zoj x o) rather than on the Gunion- 
Bertsch problem [ p2] in which the jet is replaced by a high energy beam of quarks or gluons prepared in the 
remote past. We employ as in GW [ ^, [2?j static color-screened Yukawa potentials to model interactions in 
a deconfined quark gluon plasma. The Fourier and color structure of those potentials are assumed to have 
the form 

V n = V(q n )e^ = 27r%>(q„) e -^-" T Q „ (R) ® T Q „ (n) , (4) 
where x„ is the location of the n th (heavy) target parton and 

4na s 47ra s 

where /i 2 = p 2 l± = /i 2 + q, 2 . The small transverse momentum transfer elastic cross section between the jet 
and target partons in this model is 

da el (R,T) = C R C 2 (T) \v(q)\ 2 

d 2 q d A (2tt) 2 ' 1 ' 

In our notation transverse 2D vectors are denoted as p, 3D vectors as p = (p z ,p), and four vectors by 
P = (P°,P) = [p° +p z ,p° -p z ,p\. 

The color exchange bookkeeping with the target parton n is handled by an appropriate SU(N C ) gen- 
erator, T a (n), in the d n dimensional representation of the target. (TrT a (n) = and Tr (T a (i)Tt,(J)) = 
dijS a bC2{i)di/dA- We will assume that all target partons are in the same dx dimensional representation 
with Casimir C2(T).) We denote the generators in the dn dimensional representation of the jet by a = t a 
with aa — CrI. The clastic cross section of the jet with any target parton is therefore proportional to the 
product of Casimirs, CrC2{T). 

The analytic results derived below with the reaction operator approach do not depend on the actual form 
of v, but the Yukawa form will be used for numerical estimates. Recall that in a thermally equilibrated 
medium at temperature T, the color screening mass in pQCD is given by [X = A-Ka s T 2 . Also there is a 
cut-off frequency for soft gluon modes lu p i ~ /i. We take here fj, ~ lu p i ~ 0.5 GeV for numerical estimates. In 
perturbation theory, fJ, 2 /X ~ Aira^p provides a measure of the p, the density of plasma partons weighed by 
appropriate color factors. 

This paper is organized as follows: in Sec. II we review the effect of unitarity corrections in the simple 
elastic jet scattering case. We show how the jet probability is conserved up to first order in opacity through 
the inclusion of contact double Born graphs. This is the simplest example illustrating how a factor — ^ 
arises from longitudinal momentum contour integrations in the contact limit. In Sec. Ill we proceed to the 
case of gluon bremsstrahlung associated hard probes produced in a dense medium. We discuss the basic 
diagrammatic rules and specify the assumptions and approximations used in this work. We extend the 
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algebraic classification of diagrams given in Ref. [ £7J to include virtual double Born amplitudes needed 
in the inclusive case. Sec. Ill C summarizes the rules of diagrammatic calculus that emerge from detailed 
analysis of diagrams in Appendix A through Appendix E. 

In Sec. IV the new reaction operator formalism is developed. First, operators D n , V n in Eqs. ( f76| , [Sl| ) are 
constructed from the diagrammatic rules. Products of these operators create partial sums of direct and 
virtual amplitudes from the initial hard vacuum amplitude. Those partial sums, Eq. (^), form 3™ classes 
of diagrams that can be conveniently enumerated via a tensor notation and used to construct recursion 
relations. In Sec. IV B, the reaction operator, R n = D^Dn + V n + V^, is constructed to relate the n th order 
in opacity inclusive radiation probability distribution to classes of diagrams of order n — 1. T he resulting 
simple recursion relation, Eq. (97), can be solved in closed form. The general solution, Eq. ( |101| ), is suitable 
for implementation in Monte Carlo event generators to study observable consequences of jet quenching in 
nuclear collisions. 

Color triviality of the inclusive distribution is proven to all orders algebraically with Eq. (101). The proof 
is much simpler and more transparent than in the path integral formulations [ [l^, |2^, Q and is not limited 
to quark jets. 

In Sec IV C a compact general expression for the momentum transfer averaged inclusive distributions, 
Eq. (113) is derived. Appendix F provides an independent check of this solution through second order 
starting from the amplitude iteration technique. Numerical results comparing angular distributions of gluons 
up to the first three orders in opacity are presented in Sec. V A. Analytic and numerical results for the 
angular integrated intensity distributions are compared in Sec. V B. It is shown that the induced intensity 
is dominated by the first order in opacity result that is already quadratic in L. 

A brief summary of these results up to second order in opacity was rep orted in Ref. [ [2^] . The main result 
of this paper is the new reaction operator derivation of the solutions, Eq. (101,113), that specify the inclusive 
non-abelian radiation distribution to any order in opacity. 



II. ELASTIC SCATTERING AND UNITARITY 



To illustrate how the double Born graphs cancel direct contributions to preserve unitarity we review here 
the simplest case of elastic scattering. Consider a wave packet j{p) of a parton prepared at time to and 
localized at xo = (zq,Xo) in color representation 7?. The (color matrix) amplitude to measure its momentum 
as p in the absence of final state interactions is 

M = ie ipx °j(p) x 1 . (7) 

Multiplying |Mo| 2 by the invariant one particle phase space clement e? 3 p/((27r) 3 2|p|) and taking the color 
trace gives the unperturbed inclusive distribution of jets in the wave packet: 

Consider next the effect of final state elastic interactions with an array of static potentials localized at 
x 4 = (zubi) using 

N 

Hi(t) = J d 3 x^t.(x-x 1 )T a ( l )0t(x,t)T a (i 1 , )^(i)0(x,i) , (9) 

i—l 

where D(t) — i d t and TrT a (i)Tb(j) = SijS a bC2(T)d,T / 'cLa- We will compute the three graphs in Fig. 1. The 
first order, direct amplitude to scatter with one of the (static) target partons is 

r j4 n 

M x = U?**° J jAjjiP ~ 9)A(p - q)v(q)D(2p - q) £ e^-^T a (j)T a (R) , (10) 

j=i 

where A(p) = (p 2 +ie)~ 1 and D(p) = p . The sum of double Born amplitudes in the same external potential 
is 
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M 2 = ie ipx °T a (R)T b (R) J * J* *® j(p - gi - <z 2 )A(p - qi - q 2 )D(2p - 2q 2 - qi )v( qi ) x 

N 

x A(p-q 2 )D(2p~q 2 )v(q 2 )Y / e l{qi+<!2) ^- xo) T a (j)T b (j) . (11) 




M, 



M, 




Zj Zj 



A\#\ I S q2,a 2 



FIG 1. Graphs that produce a jet with p from an initial wave packet formed at zo followed by single and 
double Born scattering center located at Z\ . 



Double Born with two different centers will not contribute because TrT a (j) = 0. To first order in opacity, 
the probability distribution of jets is given by 

d 3 N = -^Tr |M + Mi + M 2 ■ " p 



d>T ' i 2 |p|(27r)3 

= d 3 N Q + d 3 ^ + -^Tr [2 Re(Mi M*) + 2 Re(M 2 M *)] -^^TT + ' ■ ' > (12) 
«T 2|p|(27r)° 

where we separated the direct (|Mi| 2 ) and unitary correction contributions with 

ifi N f d^n d^n' 

m{2n)3 W = J (2^(2^ d nj(p-q)f(p-q')D(2p- q )D(2p-q') x 



1 i__ v(q)v*(q') CRC / T) ye^-^-'o) 

V 1 — ?V n a ^— ' 



(p-g) 2 +ie (p-q f ) 2 -ie 



(13) 



To proceed further, we consider a Yukawa potential as in Eqs. (|j|5|) and assume that all the Xj are distributed 
with the same density 

N 

p(x) = —p(z) , (14) 

where J dzp(z) = 1. Second, we assume that the observed p — (E, E, 0) = [E + , 0, 0] is high as compared to 
the potential screening scale, i.e. 

E+ ~ 2£ > ti . (15) 
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We also assume that the distance between the source and scattering centers are large compared to the 
interaction range: 

Zi ~ Zo > 1/M ■ ( 16 ) 

Finally, we assume that the source current or packet j(p) varies slowly over the range of momentum transfers 
supplied by the potential. We can then approximate the q z contour integrals ignoring the q z dependence 
of the source j(p± q). With (|l^) we can also neglect the contributions to the q z contour integrals due to 
singularities of the potentials at ±ifi±, n\ = /i 2 +q 2 . With these simplifying assumptions, only the residues 
with q z m — ie + q 2 /E + and q' z f» +ie + q' 2 /E + contribute: 



oi-i/o \3^iVi f dq z d z qdq z d z q , , < 

|P|( } ~dPp =NdR J (2tt)3 (2tt)3 (P-?)"^ 



3 -i(q-q )-(xi-x ) 



£?+g z — _ q2 _|_ j e - q/p - q' 2 — £e 

(S^ ^(P-^*^-^0( e - i(q - q ' Hxi - Xo) )«(qK(qO^^ ■ (17) 

A major simplification (that we will also use in the radiation case) occurs if the relative transverse coordi- 
nate (impact parameter) b = Xj — xo varies over a large transverse area, Aj_ , relative to the interaction area 
l//i 2 . In this case, the ensemble average over the scattering center location reduces to an impact parameter 
average as follows: 

d 2 b , s 

~a£'" (18) 

The ensemble average over the phase factor in Eq. (|l7| ) then yields 



, e 



*(q-q')- b \ — 



(2tt) 

Aj 



2 



-<5 2 (q-q') (19) 



Note that we ignored a small phase shift in Eq. ([IT]) oc exp[— q 2 (zi — zq)/E + ] w 1 on account of our high 
energy (eikonal) assumption. The ensemble averaged first order in opacity direct contribution to the jet 
distribution therefore reduces to the familiar form noting Eq. ([|) 

2|p|(27r) 3 d 3 N x N f d 2 q . . u2 da el (R,T) 
1 r|j(p - 



d 3 p A x J (2tt) 2 ' j ^ m;i d 2 q 



^\j(prjd\^ d -^^l = ^\j{p)\ 2 , (20) 

where the last line only holds if the initial packet is very wide in momentum space compared to the momentum 
transfer scale \x. 

Next we turn next to the unitary corrections. The first order term (M\Mq) — on account of Tr T a (j) = 0. 
A non-vanishing unitarity correction arises however from 

i-Tr (M 2 M*) = J ^ |^ j(p - qi - q 2 )f(p)D(2p q 2 )D(2p -q x - q 2 ) x 

V[q2) V{qi) Tr(T a (R)T b (R)) 



(p — qi - q 2 ) 2 + ie (p - q 2 ) 2 + ie 



N 1 
x{ J2e i ^ + ^- x ^—Tr(T a (j)T b {j)) ) . (21) 



j'=i 



Note that the impact parameter average constrains qi + q2 = in this case (q° = q 2 = for static 
potentials). The phase factor requires that we close the q\ z contour in the lower half-plane. We pick up 



G 



the — 2m Rcs(<7i 2 = — ie — q-2z)- This results in setting qi + q*2 = throughout. Again the residue from the 
second pole in the lower half-plane q z = —ifj,± is suppressed by the phase factor. We are left with 



1 T IM M*\ ~ ' ATrl I U M» f ^ ^ (47TQ,) 2 C R C 2 (T) 1 £+ 

— Tr (Af 2 M ) « -iNda \j(p)\ -pr-Y2~Fr\ — i j 71 — ; — 2^2 TTT 5 2~T — ' ( 22 ) 



where j_ = /i 2 + q 2 . The q 2z can be performed yielding 

dq 2z 1 £+ _ -i 1 



(23) 



faL + A? E+ <l2z - ?L - q 2 + ie 2 (q 2 + M 
The double Born "contact" contribution to the differential yield at first order in opacity is therefore 

j-Tr 2Re(M 2 M*> a -~ \j(p)\ 2 d R . (24) 

The inclusive first order in opacity elastic distribution is therefore 



cfo rf 3 p v A A ± J (27r) 2lJ ^ Myi d 2 q \X 

= \J(P)\ 2 + J dzp(z,0) J - a ei <5 2 (q)} |j(p- q)| 2 + O 0^ . (25) 

This is the first order term in the Glauber multiple scattering series. Note that when integrated over 
momentum, the double Born term cancels exactly direct contribution and therefore enforces probability 
conservation. However, for very narrow coordinate space packets, i.e. wide in momentum space so that 
\j(p~ q)| 2 ~ |i(p)| 2 ) we see that the double Born contribution actually cancels approximately the first order 
direct term differentially. 

We derived this cancelation assuming p = throughout, but it holds generally since all we needed was 
that the phase shift (zi — zo)(q — p) 2 / E + <C 1. The point of this brief review was to emphasize that the 
unitarity cancelation arose due to the factor of \ that appeared in Eq. (^) from the longitudinal momentum 
contour integration and from the factor —1 that appeared due to moving the potential from one side of a 
cut to another. This is a generic property of contact interactions that also holds in the radiation case also 
as pointed out in [ |l2| and derived here in the Appendix. 

III. SOFT GLUON RADIATION 

A. Kinematics, triple gluon vertices, and color algebra 

Consider a source J that produces a jet that subsequently radiate a gluon with four momentum k, polar- 
ization e(k), and emerges with momentum p. In light-cone components 

k= [xE + ,k~ =cu ,k] , e (fc) = [0,2^,e] , p = [(1 - x)E+,p-, p] . (26) 

Soft radiation is defined as x <C 1 so that, for example, p + ^> k + and p~ = p 2 /(l — x)E + <C k~ = ujq = 
k 2 /xE+. 

We adopt the shorthand notation of Refs. [ 

k 2 (k-q,) 2 _ (k-q t - qj ) 2 _ (k-E^q m ) 2 

Uo - 2^ ' Wi - ~^r~ ' ^ — Yuj — ' u ^ ~ 2^ ■ (27) 

In the soft eikonal kinematics, that we consider 

£+»fc+»o; (i ... j) »^P±^ . (28) 
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The soft gluon and soft rescattering approximations will allow us to simplify analytic results using: 



J( P - Q + k) 
e„(k)(p - Q + k)» 

Hp - Q) 
A{p-Q)A{ P -Q + k) 

Scattering of the jet or radiated gluon with the potential centered at position x„ introduces an integration 
/ d i q n /{2-K) i in the amplitude over the potential V(q n ). We use the GW model potential Eqs. (|,|). If 
the scattering involves a momentum exchange with the high energy jet, the vertex factor is simply —iE + in 
the eikonal limit. Neglecting the spin of the jet parton, each intermediate jet propagator brings in a factor 
+iA(p — Q), where A(p) — l/(p 2 + ie). The gluon emission vertex then gives in this approximation a factor 
+ig s (2p + k — 2Q)^, where Q is the sum of the subsequent momentum transfers. 

For scattering of the radiated gluon, consider the triple gluon tensor, r a ° 7 (fc, q) for a g(k—q, a, a)+V(q, (3 = 
0,6) — > g{k,^,c). In our model, = 0, \q nz \ = ~ \ln\ ^ lln|j an d thus the q n correspond to small 
(almost) transverse momentum transfers. If the last interaction of the gluon is with center m, then that last 
vertex is 

r Q °T(fc, q m , k - q m ) = 2ug a ~< - (k + q m ) a g°~< + (2q m - k)~< g 0a . (30) 

The color factor for the above "clockwise" kinematic convention is —f cba . 
The contraction of Eq. (^) with the final polarization gives 

T a (k; q m ) = T a ^(k, g m )e 7 (fc) = 2uje a - (k + q m ) a e° + 2g 0a (q m e) . (31) 

For the case of two gluon rescatterings the two vertices combine 

T Q (fc; q n ,q m ) = r QM (/c - q m ,q n )g^ v Y v ' (k\ q m ) 

T a (fc; q n , Qm) = 4uj 2 e a - uj(3k + q n + q m ) a e° + ff 0a (4w(g„ + q m )e) 

-g 0a (2q n (k - q m ) + q 2 m )e Q - 2(fe - q m + q n ) a {q m e) . (32) 

Note that using +g flL ' to combine the vertices above requires us to insert — iA(k — Q) for each gluon propa- 
gator, where use A(p) = l/(p 2 + ie). 

The gluon vertex factors contracted by emission vertex factors then give 

r m = (2p + k - q m ) a T a (k; q m ) « 2E+{e ■ (k - q m )) , 

r mn = (2p + k- q m - q n )c£ a {k]<ln,q m ) ~ 2E + k + (e ■ (k - q m - q„)) . (33) 

We neglected above corrections of 0(x). Eq.(|3^) also applies approximately if the emission occurs while the 
jet is off shell with p — Q, with — (0, Sq) total momentum transfer still to be acquired from subsequent 
jet re-interactions. Replacing p by p — Q adds for example 2(k + ) 2 e ■ q which is again 0(x = k + /E + ) smaller 
than the terms retained above. 

The full triple glue vertices including coupling and color algebra for producing a final color c and initial 
color a followed by color potential interactions a m and a m are then given by 

A m = T m (-r™ a ){ig s t a ){T am {m)) 

w -2g s E + e ■ (k - q m )[c, a m ]T am (m) , 

A„ m = T mn (-r- e )(-r a - a )(zg s t a )(T an (n))(T am (m)) 

^-2ig s E + k + e-(k-( il -q2)[[c,a n },a m }(Ta n {n)Ta m {m)) . (34) 

To complete the Feynman rules for our problem, we note that there is a factor 

i J{p + k - Q) e *(p+fc-Q)zo 



J(p + k)^J(p) 
ek 



kp~ ^- , 

2x 

y(A(p-Q)-A(p-Q + k)) 



(29) 



<s 



for the jet production vertex at Xq if the net subsequent target momentum transfer to the jet plus gluon 
system is Q. The hard jet radiation amplitude to emit a gluon with momentum, polarization, and color 
(k,e, c) without final state interactions is then 

M = U{p + k)e i( - p+k, > X0 (ig s )(2p + %e"(fr)*A(p + k)c 

« J(p + k)e^ +k >°{-2ig s )^c « J(p)e ipx ° {-2ig s )^- e iw °*° c , (35) 
corresponding to Eqs. (9,10) in Ref. [ p7| . 



B. Graphical shorthand 



In Ref. [ £6| we introduced the following notation, M ns!m j for the direct (tagged) n s scattering center 
amplitude for emitting a gluon between points z m and z m +\ with a final state interaction pattern encoded 
by I = J2i s &i2 t ~ m , with <7j = 0, 1 depending on whether the interaction was with the jet of the gluon. To 
include virtual double Born corrections, it convenient to associate a more specific graphical coding to keep 
track of the many different possibilities.. 

A graph consists of a gluon emission vertex, G m , for emitting a gluon between centers z m and z m+1 
as before, but also a specific set of direct interactions, A"i, CTi with <7j = 0, 1 for center i, and double Born 
interactions denoted, Oj jU , for center j. The index a,j — 0, 1, 2 denotes a contact interaction at center j with 
the jet, gluon, and both jet+gluon respectively. In this notation 



i=0 



G n 



n 



X 



(36) 



here (jj = 0, 1 is the j binary bit in the expansion of I x 2 r 



For example, the graph corresponding to 



gluon emission before the first scattering center followed by a jet interaction with center "1" and the gluon 
rescattering at center "2" is given by Mi^.i = GqXi, 0X2,1- This is simply an algebraic way of representing 
the corresponding graph. Similarly, Mo,o,o = Go and Mi o,i = Go-X^i- This notation is particularly useful to 
specify more complicated diagrams that arise from multiple direct and double Born interactions. It extends 
the algebraic approach to the topologically distinct graphs of Ref. [ by the inclusion of virtual corrections. 
Any diagram can be written in the form 



M = 



i=0 



G, 



n 



T 



3,0: 



=m+l 



T, 



(37) 



While Eq. ([37]) can be used to enumerate all "single gluon emission with rescatterings" diagrams arising 
from a target with n aligned centers, it will prove convenient to group diagrams into classes of graphs that can 
be iteratively built from a diagrammatic kernel. The two important kernels correspond to the two distinct 
possibilities for preparing an off-shell parton. In the first case the jet originates asymptotically at infinity 
and then enters the medium. It needs at least one inelastic interaction to get off-shell. This is the more 
extensively studied Gunion-Bertsch (GB) limit [ ET 
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Depending on whether the first interaction in 



the plasma is real or virtual 



D 



(R) 



Ker 



(GB) _ 



-2ig s e- (j? - (k-^V ) 



B 



(V) _ 



-2ig s f £•( 



k-qi 

k 7 (k-qi) 5 



C,Ol , 



Real initiator 
Virtual initiator 



(38) 



is the effective color current [ ^ [T^] (here given in the eikonal approximation in the small x limit). In Eq. (|38 
c and a\ are the color matrices of the radiated gluon and the soft momentum transfer. 

In the second and more experimentally relevant case, the parton is produced inside the medium (e.g., 
through A + A — * q + q + X with high Q 2 = E +2 ) with accompanying gluon radiation. The kernel or 
initiator for the hard production vertex is 
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Ker (ff) = Go = -2ig t 



ek 



(39) 



where c the color (generator) of the radiated the gluon. Classes for both Hard and Gunion-Bertsch cases 
can be constructed starting from the appropriate kernel. We focus here the hard jet case relevant to nuclear 
collisions with the hard production vertex localized at zq, as in Eq. ((39|). We consider the effect of final state 
interactions at positions z\ > zq along the direction of the jet. 

Let A denote a class of graphs with members in which the last interaction has occured at position 
Zj < Zi. Each class includes one special graph, denoted Aq, where the gluon is emitted after all real and 
virtual interactions in that class. We can enlarge this class of graphs to include a direct hit at Zi and label it 
by ADi, where Di specifies the direct insertion iteration. This new class contains 2N^ diagrams AX i and 
AXi^i obtained by a direct interaction of either the jet or the previously emitted gluon. In addition, there is 
new special diagram, AG -1 Xi^Gi = (ADi)q. where all interactions (direct or virtual) including the one at 
Zi are with the jet and the gluon is emitted after all interactions at Zi < z < oo. (G _1 amputates the gluon 
emission vertex of Aq.) The new class with 2Na + 1 graphs is constructed operationally as 



A 



AD, = AX lfi + AX it i + AG-^XiflGi 



(40) 



Starting with A = Go, ADi ■ ■ ■ D Hs is then the set of 2" s+1 — 1 direct interactions with centers "1" through 
"n s " are constructed in [ p7| . 

Similarly we can consider the possibilities that arise from inserting a double Born hit at location Z{ . This 
case differs from the above only in that now a new subclass appears because one of the legs can be attached 
to the jet line and the other one to the gluon line, i.e. AOi_2 ■ The new class AVi with 3A^ + 1 graphs 
including these virtual interactions at center i are constructed operationally as 



A 



AVi = AO ifi + AO iA + AO h2 + AG^O^G, 



(41) 



The classes and conjugate classes contributing to the gluon radiation distribution in the opacity expansion 
at order (L/X) n cx (a e i/A±) n for n = 0,1,2 are listed in the table below. Each class is obtained through 
ordered insertions from the hard kernel Ker' ff ' indicated by Go- 

1. (L/X)° oc (a e i/A±)° contributions 

• 1 (n s = 0) x (n s = 0)t 



GqGJ 1 graph 



(42) 



2. (L/X) 1 cx (cr e i/A±) 1 contributions 

• 4 (n s = 2) x (n s = 0)t 

• 9 {n s = 1) x {n s = l)t 



Go^lG 4 graphs 



GoD^Gl 9 graphs 



(43) 



(44) 



3. (L/X) 2 oc (a e i/A±) 2 contributions 

• 13 (n s = 4) x (n s = 0) f 

• 57 (n s = 3) x (n s = l) 1 



G ViV 2 Gl 13 graphs 



GoD^DlGl 30 graphs 
G q ViD 2 dIgI 27 graphs 



(45) 



(46) 
(47) 
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• 65 (n s = 2) X (n s = 2)t 



G DxD 2 dIdIgI 49 graphs (48) 
GoVxV^Gl 16 graphs (49) 

The left column sorts contributions in terms of diagrams with fixed number of interactions n s ,n' s in the 
amplitude and complex conjugate amplitude. For example, there are 57 contributions involving amplitudes 
with 3 interactions times complex conjugate aplitudes with 1 interaction that contribute to second order in 
opacity. Since two of the three interactions in the amplitude must involve the same center, these 57 diagrams 
subdivide naturally into the two class structures (^6|,^) indicated on the right column. We note that (fl2"|), 
( fi"4| ) and (^|) correspond to the direct diagrams computed in [ |27j for the exclusive tagged spectrum. 



C. Diagrammatic calculus 

Whereas the time-ordered perturbation techniques [ [l(J used in [ |2lJ greatly facilitated the listing of higher 
order contributions, it is useful to show how those shortcuts emerged from Feynman diagrams especially to 
check the shortcut rules that apply to the double Born contributions. Detailed diagramatic evaluation from 
which those rules emerge are presented in Appendix A through Appendix E. They exhaust all possibilities 
that we can encounter for single and double Born interactions in a potential model with multiple centers. It 
is important to note that a double Born interaction at position Zi can be thought of as the "contact" limit 
of two direct hits at positions z\ and Zi+i respectively 

AOiai oc lim AXi, bt X i+ i, b . (50) 

We summarize here the kinematical and color parts of this diagrammatic calculus: 

• In all diagrams the jet production factors out in the form 

J(p) e tpx ° 

under the soft-gluon soft-interaction assumption Eq. (p9j). Therefore the jet distribution can be factored 
out under the assumption that J is broad. 

• For graphs involving direct interactions only, i.e. of class GqD\D^ ■•■£){•■■, the Feynman diagrams 
reproduce the time ordered perturbation results listed in Ref. [ ^7j. The simplest one interaction 
example is discussed in Appendix A. We note that although Appendix B through Appendix E are aimed 
at clarifying the contact limits, their intermediate results in the well-separated zi+i — Zi = X 1/fi 



case are directly comparable to [ 27 



• A diagram which involves a contact contribution of the form AOi$ or AO^x, i.e. both interactions are 
on a single line, is kinematically identical to a diagram in the well-separated case , oc XifiXi+x,o and 



oc X^iJQ+i,! respectively, but with Zj = Zj+x (50). This can be seen in from Appendix B for the final 
state virtual interaction of the radiated gluon and in Appendix C for a virtual final state scattering 
of the jet. There is, however, a multiplicative factor of | per contact interaction that emerges from 



contour integrals. Eq. (B8) shows how the strength interpolates between 1 and | as the separation 



varies for a Yukawa potential. Eq. (B15) shows the generality of the factor | in the contact limit. 



• For the case when one of the legs in the contact contribution is attached to the jet line and the other 
to the gluon line, i.e. .40^2, the amplitude is kinematically identical to a diag ram AX^qX^+x^x in 
the well-separated case with z% = Zj+i (|50|). There are no additional numerical factors for such mixed 
contact term as shown in Appendix D. We also note that the contact limit of the diagram oc Xi^Xi+x.o 
should not be added to avoid over counting (i.e, they are topologically identical in the contact limit). 

• We show in Appendix E that the class of diagrams that have the gluon emission vertex between the 
two momentum transfers of a double Born interaction are suppressed in the "contact" limit by factors 
0(k + 1 E + ) and 0(/j,/k + ) and thus can be neglected in the framework of our approximations (p8|). 
This verifies the naive argument from time-ordered perturbation theory that such diagrams are of 
zero- measure and vanish in the contact limit. 
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• Each power of the potential comes with a weight (elastic scattering amplitude) 

H) / ^ V (0, qi )e-^ b -T ai (z) , 

that factors out as well under the assumptions ( |29| , p8| ) and is still to be integrated over. The amplitudes 
above eventually build into the elastic scattering cross section oc a e i . 

• The appropriate color factor for the jet+gluon part can be constructed following Ref. [ . In the case 
of a virtual interaction certain simplifications, namely 

a-iO-i = C R 1, ai[c,a,i] = — ^-c , (51) 

can occur even at the amplitude level. 

• The last critical rule that emerges is that direct terms come with a plus sign while each virtual interaction 
contributes a minus. 

These rules are illustrated to compute the classes of amplitudes needed up to second order in opacity in 
Appendix F. 



D. Impact parameter averages 

Consider next the averaging over the transverse impact parameters = x,; — xo. A typical term that 
appears to second order in opacity is GqXi O2^X[ G . This involves one direct interaction with center 
"1" and a double Born with center "2". We consider separately the average over involving direct and 
virtual interaction centers. 

For a scattering center, i, that appears in a direct interaction the impact parameter average takes the form 

<-W = (-/ ^H) / §%v(0,<u)e-^(+i) J **v* M )e^< ) 
d 2 q, d 2 q> (2^ 2 £ 2 ( qi -q^ 



-«(0,qi)«(0,q<) 



(2tt) 2 (2tt) 2 A± 

= • • • ihS wF Hcii)? I d2q ^ 2(q,; _ q5) " ' ' (52) 

where \v\ 2 is defined as the normalized distribution of momentum transfers from the scattering centers. In 
terms of Eq. (||), it is given by 

a e i d 2 q 7T (q 2 + fl 2 ) 2 
where in the Yukawa example, the normalization depends on the kinematic bounds through 

111 

(54) 



Me// M 2 Qmax + M 2 

and insures that f qmax d 2 q|i)(q)| 2 = 1. In numerical estimates we take q 2 la2 . = s/4 w 3Efi. For asymptotic 
energies /^ e // ~ f 1 - Note also that unlike v, the barred v is independent of color factors and thus applies to 
both quark and gluon rescattering. 

If the interaction with center i is a double Born rather than a direct one, then the only change relative 
to ( j52] ) is that both q^ and qj come with the same (±i) 2 = —1 rather than the (—«)(+*) = 1 factor and the 
phase shifts change relative sign. If the double Born interaction is in the amplitude, as in the example noted 
above, then the impact parameter average leads to 
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■)a ± 



d z hi 



H) 



cPq, 



= •••(-1) 



A i v "'J (2tt) 2 

d 2 q, d 2 q^ (27r) 2 <5 2 (q 4 + q' J ) 



H) 



(2tt) 2 (2tt) 2 
r2. 



-4 



(27T) 2 

«(0,qO«(0,q<) 



«(0,q-)e 



-iq/ hi 



(-l)^/^l-(q 4 )| 2 / d 2 q^ 2 (q l +q' i ) 



(55) 



The same holds true if the double Born is in the conjugate amplitude. The results above is the basis for the 
last graphical rule of the last section. A contribution to n th opacity involving m "virtual" hits produces a 
factor (-l) m (a/A±) n . 

A key to the analytic derivation in the next section is that the transverse impact parameter averages in 
Eqs. j5^j55| ) diagonalize the products of amplitudes and complex conjugate amplitudes in the q^, qj variables. 

The question of which cross section a enters above depends on the details of the color algebra. At this point 
it is not at all clear since any particular contribution to a given order in opacity has an extremely complicated 
color structure in general [|27J]. However, as we prove in the next section the color algebra simplifies drastically 
in the sum over all possible real and virtual contributions [ [ll], . The color "triviality" of single inclusive 
distributions proved in the next section via the reaction operator approach, states that after summing all 
contributions, each integration over q^ is accompanied simply by a factor ( C^Ci (T)/(1a)- This implies that 
a = <j g is the gluon elastic jet cross section that appears in both Eqs. (|5^ , |55| ), rather than the jet one 
(cx C R C 2 {T)) . Recall that 



Ca 
Cr 



<Tel(R,T) 



and hence 



C^An = CrX 



(56) 



(57) 



Color triviality implies that the opacity expansion is in terms of powers of the gluon rather than the jet 
elastic cross section, i.e. oc (L/X g ) n . 



IV. RECURSIVE METHOD FOR OPACITY EXPANSION 



The inclusive gluon distribution to ©((L/A) 1 cx (<J e i/A±) 1 ) is a sum of 3 2 direct and 2x4 virtual cut 
diagrams. If we want to proceed to second order in opacity 0((L/X) 2 oc (a e i/A±) 2 ) there are 7 2 direct plus 
2 x 86 contributions. It is therefore useful to develop an recursive procedure for writing the sum of amplitudes 



in a certain class of diagrams. We will ignore the jet and elastic scatterin 
us with simpler amplitudes similar to the ones discussed in Refs. [ |26|, f!7 
Hard, Gunion-Bertsch and Cascade terms 



aarts of the diagrams which leaves 
We recall the definitions of the 



k 

k 2 ' 



H 

B; = H — Cj , 



(k q„ - q l: 



B(i 1 i 2 -"£m)(il.J2 ■■■in) ^ 



(k - q n - q l2 - 

3n. 



-q*J 2 



(58) 



A. Amplitude iteration 



There are two basic iteration steps that one has to consider to construct the inclusive distribution of gluon. 
The first one represents the addition of a "direct" interaction D n that changes the color or momentum of 
the target parton located at illustrated by Fig. 2. 
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+ e 1 



(COo-Mi) Z; 




Elj Col(c) 




k — > k - q ; 



Col(c -» [c,a;]) 



- (-^Bie^^tcaJTd 



FIG 2. Diagrammatic representation of the sum of amplitudes generated by ADi = AXi,» + AXi,i + 
AG 1 Xi,oGi. 



The second iteration step corresponds to a double Born "virtual" interaction V n that leaves both the color 
and momentum of the target parton unchanged and is illustrated in Fig. 3. 

Our goal here is to construct operators D n , V n that can be used to construct recursively the new amplitudes. 
It is important to note that each new diagram produces a difference of two phase factors with the exception 
of one special diagram, (AD„)q = AG~ 1 X„ t oG n , corresponding to the emission of the gluon after, z n . The 
upper limits of emission is z — oo in that case, and only the lower bound phase factor, — e luJ ° z ™ survives 
on account of the adiabatic damping factor exp(— e\z\). Therefore all but one diagram can be divided into 
two parts which can then be recombined into more easily interpretable Hard, Gunion-Bertsch and Cascade 

"271 - 



amplitudes (pq) as in [ 26 



The sum of amplitudes in class A can be denoted by 



A(x, k, c) = y~] Aa(x, h)Col(c) a 



(59) 



where A a (x,Js) represents the kinematical part, Col(c) a stands for the color matrix for the distinct graphs 
in this class enumerated by a. There is of course considerable freedom in rearranging the terms of this 
sum. Since by definition classes are constructed by repeated operations of either of one of three operations, 
i,Di,Vi it is convenient to enumerate the 3™ different classes of diagrams via a tensor notation, Ai t ...i n , 
where the indices ij = 0, 1, 2 specify whether there was no, a direct, or a virtual interaction with the target 
parton at Zj. The sum of amplitudes for class, -4j 1 ...j„, can be written explicitly in the form 



A 1 ...»„(x,k,c) = Y[ 0>o,i m + Si,i m D m + 5 2: t m V m ) G (x,k, c) 

111 — 1 



(60) 



Here Go is the color matrix amplitude that corresponds to the initial source or kernel of the jet and gluon 
in the limit of no final state interactions. 

The inclusive induced "probability" distribution at order n in the opacity expansion can then be computed 
from the following sum of products over the 3" classes that contribute at that order: 



P n {x,k) =A it - in (c)A il ... in (c) =Tr^2 ■•• A 1 ..^,k 1 c)A il ... in (x,k,c) , 

i 1= i„=0 

where the unique complementary class that contracts with Ai t ...i n is defined by 



(61) 
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Ai 1 ...i n (x,k.,c)= Y\ \So,i m %i + Si iim D m + S2,i m ]Go(x,'k,c) , 

ni—1 

n 

A n -^(x,k,c) = Gl(x,k : c) JT (s , lm Vl +Si,i m Dl + d 2< i m ) • (62) 

m— 1 

For example, .A0210 = D3V2G0 for ^2012 = V^DsViGq. Note that AA includes products of classes and 
complementary classes with different powers of the external coupling, g s , but every product in the sum is 
the same order, a 2 s n+1 . These mixed terms are the unitarity corrections to the direct Ar"' 1 A\...\ term that 
contribute to inclusive processes when the target recoils are not tagged, i.e. measured exclusively. Note that 
P n is not positive definite except for n = 0, and we refer to it as a "probability" only figuratively. 

With this tensor classification and construction, it becomes possible to construct P n recursively from lower 
rank (opacity) classes through the insertion of a "reaction" operator as follows: 

P n = A il " 4n - 1 R n A il ...i n _ 1 , Rn = Dlfin + Vn + ■ (63) 

We emphasize that it is possible to write the inclusive probability in this simple form only because the trans- 
verse impact parameter averages in Eqs. ( |52"| , [55"| ) diagonalize the products of amplitudes and complementary 
conjugate amplitudes in the q^, variables. If the transverse area of the target were not large in comparison 
to the cross section, then off diagonal components in those variables would appear that would have to folded 
with a suitable transverse Glauber profile functions as discussed in [ [|. The simple algebraic structure of 
the problem in Eq. ( |63| ) therefore hinges on the validity of the transverse averaging via Eqs. ( 

To construct D n , consider first the contribution AX n fl from Eq. (flo[). In this case, a new direct interaction 
is attached to the jet line at z n . For the partial sum of amplitudes in Ai 1 ...i n _ 1 corresponding to emitting the 
gluon before the last real or virtual interaction at Zf < z n _i, the direct interaction of the jet at z n simply 
multiplies that partial sum by a n . No extra phases are introduced since the energy of the jet is assumed to 
be very high. As noted before, for every class, there is a term in Ai 1 ...i n _ 1 corresponding to a graph where 
the gluon is emitted after the last interaction point Zf. That term has the explicit form 

(-i)^ ,( ^ 1 - , "- l) H(- e ^)cT e/ (A 1 ... l _ 1 ) , (64) 

where the elastic color matrices for that class and its complementary are denoted by 

TeiiAi,...^) = (a^i) 1 - 1 • • • M 11 , rJjCtf 1 -*"- 1 ) = (ai) 2 - 11 • • • (an-i) 2 - 1 - 1 . (65) 

For example, T e i(Aio2i = D^V^DiGq) = 04(0303)01 = Cna^cLi. Eq. ( ]65| ) implies the following identity holds 
for elastic color factors: 

Tl(A i ^ i ^)T el (A il ... in _ 1 ) = C^ 1 (66) 

valid for all 3™ _1 components of rank n — 1 classes. 

The factor {—^) Nv in Eq. ( |64| ) arises because every virtual contact interaction in the class acquires a factor 
— i from the the contact limit of the contour integration over longitudinal momentum (see Sec. II and the 
Appendix). The numbers, N V ,N V , of such contact interactions in class, Ai 1 ...i n _ 1 , and its complementary 
class, A 11 '" 1 "- 1 , are 

n— 1 n— 1 

m—1 m— 1 

For example N v = 2,N V = 2 for ^01202 = V5V3D2G0. A useful bookkeeping identity that we will need is 

E H)^"- 3 (-\) N ^- ] = H - \ + i) m = • (68) 

Now consider how a jet interaction at z n changes the graph Eq. (|64[). Besides being multiplied by a n , the 
only other change is that the phase factor changes into a difference of phase factors as 



52 H). 
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The part, oc e luJ ° Zf , on the right hand side can be added back to the sum of all the other graphs with 
emission before Zf. The modified sum a n Ai 1 ...i n _ 1 therefore includes this part of Eq. ([39|). In summary, the 
component of the D n operator that takes into account the jet rescattering at z n is simply multiplication by 
the color matrix, a n , in the jet representation. This then identifies the first part of the D n operator as 

D n = la n + --- . (70) 

The left over part oc e 4 "° z " of Eq. ( |69| ) can be combined with the graph (AG~ 1 X rlj0 G n ), e.g., Mi^o, M2,2,o 
in Appendix A and C, that accounts for gluon emission after z n . This combination leads to the following 
extra term in A^-.-i^: 

.(.l^.^-^) Hgia**, [ca^T^...^) . (71) 

The sum of graphs A{X nfi + G^ 1 X rhQ G n ) is therefore a n Ai 1 ...i n plus Eq. (|7l|). 

To finish constructing D n , we need to consider the effect of an elastic scattering of the gluon produced at 
an earlier point in the medium. The gluon "cascade" interaction graph AX nj \ (see e.g., A/2,0,3 i n Appendix 
B) introduces the following changes in the amplitude: First an extra phase shift exp[i(uJo — aj n )z n ] arises due 
to the difference of gluon energies before and after the transverse momentum exchange, q„, at z n . This is in 
contrast to the jet rescattering where that energy change was negligible. Second, the gluon color is rotated 
by c — > [c, a n ]. Finally, the transverse momentum in A(x, k, c) is shifted backward by k — > k — q„ in order 
that same final state k is reached. This part of the D n operator therefore simply transforms 

Av-i„-i(M S n A h ... lrl _ 1 = e^o-^nK A h ... lrl _ 1 (k- q n ,[c,a n }) . (72) 
The "shift" or gluon scattering operator that implements gluon rescattering at z n is thus defined by 

S n = if cand x e <(<"o-«.0*» e <q»-6 _ (73) 

Here, b = i Vk is the impact parameter operator conjugate to k such that 

e^ 6 /(k)-/(k-q) . 

This transverse momentum shift together with a class independent phase shift and a color rotation (via 
the structure constants f abc ) enlarges Ai 1 ...i n _ 1 to include all but one of the gluon final state interaction 
diagrams with an interaction at z n . 

The left over graph not included in (j7^) is the one where the jet scatters elastically up to previous last 
interaction point Zf, emits a gluon with transverse momentum k — q„ between Zf and z n , and the gluon 
scatters at z n (see e.g., Mi, 0,1 in Appendix A). The integral over the emission interval leaves a characteristic 
differences of phases e lUnZn — e luJnZf , while the elastic scattering introduces a phase shift factor e l (^o-"»i)z„ 
as well. The kinematic part of this cascade amplitude is C„. The part proportional to e luJnZn contributes 
therefore the following extra term 

+ (-i) A '"C„ e ^-" 1 ) 2 "e^ z "[c,a„]T e/ (A 1 ... l „_ 1 ) , (74) 

where again N v is the number of contact interactions before z n . Recalling that B„ = H — C n has the 
physical interpretation of a Gunion-Bertsch gluon source term at z n , the two extra terms from Eq. (j74|) and 
Eq. ( fn] ) can be conveniently combined to form 

BnAir-i^M = - (-\) N ^ -^\r- d x B„e"» z » dT^A^..^) 

= - n (- 2 J B - e " oz " [c ' an] a * zi ■ • ■ ^ ■ (75) 

Combination of three steps X ni o+X nt x+G X n fiG can therefore be summarized by the following operation 
on rank n — 1 class elements: 
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■Dn^-ii ■■■i T1 _i(^-;k, c) — (o-n H~ S n -\- B n ■■■i n -\ k, c) 



-1) 



B„ e v 



[c, a„]T e j (Ai- 



(76) 



We note that the common factor to all diagrams 2ig s e •(•••) is suppressed for clarity. 

We now turn to the second class of diagrams, V n Ai 1 ...i n _ 1 created by inserting contact interaction at z n . 
This new class is constructed via Eq. ( [ill) an( A illustrated in Fig. 3. 




Col(c) I 1 




Cr+C a 



Col(c) 



^(coo-cO;) zj 




a ; Col(c — > [c,aj) 



-7 (-7)^ 



e c T, 



cl 



FIG 3. Diagrammatic representation of the sum of amplitudes, generated by AVi 
AOi^+AG^O.nG,. 



AO lfi + AOiA + 



For a double Born interaction at z n , the transverse coordinate integration can be performed before multiply- 
ing by a complementary hermitian conjugate amplitudes since the coordinates of a virtual interaction cannot 



occur in the complementary class. The transverse coordinate integration on the center of the i th interaction 
enforces (see Eq. (|l9|)) that the net momentum transfered in the two legs vanishes via a factor <5 2 (qi + qj)- 
Consider first the sum of graphs where the two legs of the contact interaction are attached to the jet line 
(AO n ,o) (see, e.g. M| 00 in Appendix C). The situation is similar to Eq. (70) except that two color factors 
o-ndn = Cr multiply the rank n — 1 amplitudes instead of one. In addition a factor of — i appears due to 
the contact limit of the contour integration. Thus, this part of the virtual interaction simply transforms 



Ai 



o„ 



Cr 



(77) 



The second term, AO n ^\ in Eq. ( |4l| ) is the contact with the gluon at z n (see M$ 3 in Appendix B). Since 
^ca n d^jda n e _ Cjjj the color structure of this term is also trivial. In addition, since the total momentum 
transfer is now zero, no extra phase factor arises and the gluon momentum in Ai 1 ..-i n _ 1 remains unchanged. 
The contact gluon interaction at z n therefore contributes 



Ai 1 ■■■i n _ 1 



0„,i 



C A 



•Ail ■ 



(78) 



The mixed contact interaction AO 1.2 (see, e.g. M| 01 in Appendix D) is the most subtle. Even though 
the net momentum transfer vanishes, in this case both the jet and gluons receive (opposite) non-vanishing 
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transverse momentum transfers, ±q„. This leads to a highly non-local modification of . The contact 

limit of the longitudinal momentum transfer contour integration produces in this case a factor of —1 as shown 
in the Appendix. The elastic scattering of the gluon produces multiplicative phase factor, e. l ^ J °~' JJ ' n ' Zn . The 
jet elastic scattering introduces no new phase in the small x limit as before. The transverse momentum 
exchange, to the gluon requires that k is again shifted back by q n as in Eq. (|7^). The gluon color is also 
rotated as in (|7^). Furthermore, the scattering of the jet introduces another color matrix a n . This mixed 
contact interaction therefore simply produces the following variant of (J7^): 

Aii---i n — i (k, c) -> a n S n A{ 1 . . -i n _ i — & 1 nAi 1 ■••i n _i (k q n ,[c,tt n ]) , (^9) 

with again one extra term that can be recombined with the amplitude of gluon radiation after the contact 
interaction at z n to form a virtual Gunion-Bertsch source 

- a n B n A il ... in _ 1 = - (-^^-^ ^B n e^cT el (A n ... ln _ 1 ) , (80) 

where we used a n [c, a n ] = — ^Cac to simplify the color algebra. 

The double Born interaction at z n can therefore be implemented by the following operator 

V n = —\{Ca + Cr) — a n S n — a n B n = —a n D n — \{Ca — Cr) . (81) 

Specifically, the virtual iteration of a rank n — 1 class gives 

^nAi- l -"i n _i k, c) — — x^ii •■ - z n - 1 (*£; k, c) 6 ^ ^ (2 nw A.^ 1 ...^ n _ 1 (x, k q n , [c, ) 

-(-i) JV "%B n e~cat- 1 I ---ai 1 . (82) 



B. Reaction operator recursion to all orders in opacity 



We emphasize that this recursive process is a completely general and applies both the hard emission kernel 
Ker (H) Eq. (||) and the asymptotic Gunion-Bertsch Ker( GB ) Eq. (||). 

The key to unraveling the interplay between direct and virtual interaction is the operator identity Eq. ( |8l] ) 
between D n and V n - This makes it possible to relate the n th order in opacity gluon emission "probability" 
distribution P n to the probability at (n — l) th order by expressing the reaction operator R n in Eq. (|63|) as 



Rn = (D n - a n )\D n - a n ) -G A = (S n + B n )\S n + B n ) - C A ■ (83) 

A further major simplification occurs because both S and B involve the same gluon color rotation through 
if cand . We show next how this unravels the color algebra and reduces it to simply powers of Ca and Cr. We 
therefore prove color "triviality" of the reaction operator, a property that is implicit in the work of Rc fs. [ 
[TT| , |l3| and which has been proved in a somewhat more involved way via path integral techniques in [ 34 
for quark jets only. 

First, we note that the shift operator has the property that 

SlS n = Oie- q " Vk e- q " Vk . (84) 

Therefore S^S is diagonal in color space. This allows us to compute the effect of an extra gluon cascade 
interactions through simple transverse momentum shifts 

A^-^-'iStSn - Ca)A 1x . ..<„_!= C A (^"^^(k - q„, c)A n . (k - q„, c) - ^ 1 - 4 - 1 (k,c)A 1 ... 4 „_ 1 (k,c)) 

= C A (Pn-i(k- q„) - P„-i(k)) = C A (e iq ™' 6 - l) P„-i(k) . (85) 

Next, consider the modulus square of the extra Gunion-Bertsch source for two or more scatterings of the 
jet+gluon system (n > 1): 
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A^-^BtBnAi^i^ = C A |B n | 2 J2 H) ("I) ^(^ 1 "'^ 1 )ccT e; (A 1 ...„_ 1 ) 



ti,--,»n-i 



n-1 



caCr ib„i 2 x: n (-i) d2 - im K) Vim w 51 " 



(86) 



21 .■■■,i Tl — 1 m—1 



where we used Eq. (|66|), and the bookkeeping identity (68). Therefore, the virtual corrections cancel the 
direct extra Gunion-Bertsch source term for n > 1. For the special case with n = 1, however, this diagonal 
term survives 



Gj J B 1 t B 1 Go=C A C R |B 1 | 2 . 

Finally, we have to calculate the interference between the cascade and Gunion-Bertsch source terms in ([i3]) . 
The two gluon color factor f ca ™ d in Eqs. (73 7^) again contract to form a color diagonal factor Ca- This 
leaves 



(87) 



2YLeA ll -^- 1 B} l S n A ll ... ln _ 1 = -2C A B n • (Re e 
where the I n _i term in the iteration step is easily seen to be 

In-l = ( — 5) T^,[(A 1 " 1 )cAi 1 ..-i n _ 1 



E {(-5) 5q ' i "- i ^ 1 -;"} c (v^-i +^i^- 1 a,-i + ^„_ 1 k-i)a 1 ...^ 



in -1=0 



E (-i)^" 4 ' 1 " 1 " 2) r e t i (Ai---i„_ 2 )(-5 a «-i a n-iC + a„_ 1 c J D„_ 1 + cV r „_i).A il . 



(89) 



This contribution builds up on the class of diagrams Ai 1 ...i n _ 2 under the action of an iteration operator 
which can be reduced using Eqs. ([73 81) into the form 

- \a n a n c + a n cD n + cV n = [a n , c]D n — c = [a n , c](S n + B n ) - Cac 



CA ( e i(u,o-u n )z ne i* n -b _ i) c+ [ aniC ] Bn 



(90) 



Next, we note that the contributions from B actually cancel exactly for n > 2 in the sum over all classes 
due to the bookkeeping identity: 



2) ^eliA 1 n 2 ) [ a n — 1 7 c\B n — \Ai-t ...i n _ 2 



E 

.■■■in-2 

-CUC fl B n _ie <w °*>- 



1 E 



l^iV„(«4 il ... i „_ 2 ) ( _ i\N v (A il ... in _ 2 ) 



(-|) JVo(Al '" < "- 2) ^(A 1 ... in _ 2 )T ei (A 1 ...i„_ 2 )=0 . (91) 



For the special case with n = 2 we can evaluate Ii directly from the definition Eq. ( jS9|) and the initial hard 
amplitude to obtain 



Ii = C A (e 



l)cG + [01, c]BxG = -C A C R (( e ^o-^i)* V qi-b _ ^h^o*, + Bie i 



(92) 
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In the last line the Gunion-Bertsch amplitude was rewritten with the help of the shift operator Eq. ( |73| ) as 
being derived from the hard vertex kernel 

Bie «^o*i = _ ^giK-wO^igiqi-fi _ A He 4WoZl . (93) 

Eqs. ( |89| , p0| , ^ ) therefore imply that I„ obeys the recursion relation 

I„ = c A (e^~^ z -e^ - l)I n _! - S nA C A C R B ie 1 ^ . (94) 
With I = -C R He iuJ ° z ° and Eq. @ for n > 1 we solve @ in a closed form 

I„ = C fl C3 1 J| ( e H^o-^ m )^ e iq m -b _ ^ I H (e iW0Zl - e M02 °) , (95) 

lm=l J 

where the product is understood as ordered for left to right in decreasing order in operators labeled by m. 

Combining Eqs. ( ^5||8^j87| , |88| , |95| ) that specify how the reaction operator acts between the contraction of 
rank n — 1 classes, wc find that the inclusive radiation probability, Eq. (|63|), obeys a simple recursion relation: 

P n (k) = C A (P»-i(k - q„) - P„_i(k)) - 2C A B„ • (Re e- l ^e^H n -i) + S^CaCr^ | 2 . (96) 

We can solve Eq. ( p6| ) with the initial condition 

P = C R H 2 . (97) 

We introduce the following notation for the separation of the scattering centers Az„ = z n — z n -\. For n = 1 
the solution is 

Pi = C A C R (C? - H 2 + B 2 + 2Bi ■ Ci cos(wiA«i)) = -2CUCRB;! ■ Ci(l - cos(wiAzi)) . (98) 
For n = 2 we need 

Re e" ia;2Z2 7i(k - q 2 ) = -C A C R Re e ~ M222 (c (12) (e l( " 2Zl ^< 12 ) A2l) - e luJ2Zl ) - C 2 (e^ 22 ° - e^ 221 )) 

= C A C R [-C (12 )(cos(o;2Az2 +u il2 )Az 1 ) - cos(w 2 Az 2 )) 

+C 2 (cos(w 2 (Azi + Az 2 ) - cos(w 2 Az 2 )) ] . (99) 
Therefore, the n = 2 distribution is given by 

P 2 = 2C\C R [ Bi • Ci(l - cos(wiA*i)) - B 2(12) • C (12) (l - cos(w( 12 )A«i)) 
+B 2 • C(i 2 )(cos(w 2 Az 2 + w (12 )A2i) - cos(w 2 Az 2 )) 

-B 2 -C 2 (cos(w 2 (Az 1 + Az 2 )-cos(w 2 Az 2 ))] . (100) 
For n > 1, we can use (^H|) to obtain the general solution for the gluon probability at n th order in opacity 

P»(k) = C A (e^ - l)P„_i(k) 

-2C fl C2B„ • Re e -^»^e iq °' 6 | ( e *(wo-"™)*» e iq '»- 6 - 1)| H(e lu » 21 - e Mo2 °) 

n j n ] 

= -1C R C\ Re ^ < JJ :<"' - 1) I Bi - e *»-6 e -*"o*« x 
i=i J 

j JJ ( e <(«o-wm)*m e iqm-b _ j)! H(e^ 0Zl - e^ 020 ) , (101) 
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where the last applied to n = 1 as w ell. This is therefore the complete solution to the problem and is the 



central result of this paper. Eq. (101) provides not only an algebraic proof of color triviality of the inclusive 
gluon distribution [P n oc CrC 1 }), but in fact gives the complete angular dependence for arbitrary Zi and q„. 
It also applies to both quark and gluon jets. The derivation through the reaction operator and recursion 
relations is furthermore much simpler and transparent than previous derivations [ [H], |l3], |34| . The result is 
also more general and versatile for applications to nuclear collisions. 

In this form, it is for example possible to im ple ment P n directly as an "after-burner" in a Monte Carlo 
event generator, e.g. ZPC [|| || or MPC [||,||. Such programs solve relativistic 3+1D transport theory 
and provide much more reali stic d istributions of target parton coordinates Zi along the path of a high energy 



jet in nuclear collisions. Eq. (101) also makes it possible to consider evolving potentials since the distribution 
of the q„ remain unspecified. For example, the screening scale that characterizes the mean q„ transfered 
at z n , could itself depend on n because of the variation of the local density at that point. In this way, jet 
quenching in realistic, evolving nuclear geometries can be calculated more accurately. 

C. Ensemble average over momentum transfers 

We now consider the integrated gluon distribution over the q m via the normalized squared scattering 
potential Eq. ( |53| ) 

„ n 

(P„(k))„ = / \{ {dq m v 2 rn (q m )} P n (k; qi .--q n ) , (102) 

m=l 

where we allow for the possibility that the effective potential could change along the path of the jet. We 
hold the z rn fixed as yet. ^From Eq. (Q) we need to consider for example 

J dq n ^(q„) (P„_i(k - q n ) - P„_i(k)) = / dq n (v 2 n (q n ) - <5 2 (q„)) P n _i(k - q n ) . (103) 

Since <5 2 (q„)B„ = 0, we can also replace w^(q n ) by v^(q n ) — S 2 (q n ) in the average over the inhomogeneous 
term. However, due to the special form of the operators defining I„ in Eq. ( |95| ) we can also rewrite the 
integral over q^ for i = 1, ■ • • , n — lin analogous form 

/n n—1 
i J {dqj v](qj)} B„ e ^^ e ^-b JJ ^i(uj Q -a, m )z m & i^ m -h _ ^ 
j=l m=\ 

r, n n 

= / fi{*fc [^(q,)-^ 2 (q,)]} B n e-^- J] e^-^-e^ , (104) 
j=l m=l 

where the product of operators is again path ordered. This product of operators can be simplified as follows 

n 

& -iu z n J~J ^(aJo-WjOzjgiqj-b _ £ -iu n z n e i(uj n -uj (n _ 1 „ } )z„_ i _ _ _ gi(w (2 . . . n) . , n) )zi giQb 
J = l 

= exp ( -iJ2u u ... n) ( Zj - z s -i) \ e -*-P--)*o e <Q-fi. , (i 05 ) 



where Q = ^2 is the total momentum transfer, and 

(k - q m q„) 2 



■'(m,-- -,n) 



2xE 



(106) 



Therefore, we can evaluate 



e -ft-o*„ Yl e <(«o-<"i)*i e *'» ,fi H ( e ^ 3i _ e iu>o*o) = C (1 ... n) e i *"'"(e ia '( 1 -"> (zi - zo) - 1) , (107) 
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where is the gluon elastic scattering phase shift from zq to z n . The partial (eikonal) phase shift 

{J z " 1 dz {k z (z) — xE)) due to gluon rescattering from z m _i to z n is given by 

m rn 

$n,m = ~ y^U(k-n)i z k ~ Zfc-l) = ~ ^ U(k-n) AZfc . (108) 
fc=l fc=l 

Eq. (107) makes it possible to write the recursion relation for the momentum transfer averaged probability 
for n > 1 in the form 

(Pn) v =C A J dq n (t£(q n ) - <5 2 (q„)) ( P„_ x (k - q„) >„ 

/n 
(^( qi )-^ 2 (q 4 ))} B„.C (1 ... n) Ree l *"-(e^-">( Zl -^-l) . (109) 
i=l 

With Pi given by Eq. (|98|), the solution for n > is 

« n 

( P„ >„ = - 2C R C% J J] {dq, (« 2 (q t ) - <5 2 (q,))} 

i—1 

n 

^ B (m+1 ,... >n)(ro ... n) • C (1 ... n) Re (e**«.« (e*-ti-») (*i-*> - 1) ) . (110) 



X 

m— 1 



where for m = n, B ( „ +1 ...„ )(n) = B (0)(ri) = B (n) . 

Recall that we cannot take the contact limit z\ = zq, which apparently vanishes, because our derivation 
assumed that z m — z m -\ were larger than the range of the force 1/fi. However, it is clear from Eq. (110), 
that it is the first step phase difference 



(!...„) (Xl-Zo) _ I 

that controls the bulk of the LPM destructive interference effect that suppresses radiation when the formation 
length \/<jj(i... n \ is long compared to (z\ — zq). The subsequent interactions merely modulate this effect with 
an elastic scattering phase shift. 

Finally, we can restore the proportionality constants by multiplying with a s /-K for the production vertex, 
1/tt for the d 2 k measure, and Y[j( a gU)/Ax) along the path to convert the Vj back into Vj from (|53|). In 
addition, we must multiply by the combinatorial factor, 

AM N n , s 

(112) 



n\(N — n)! n! 

that counts the number of ways n target partons out of can be within the interaction range of the jet+gluon 
system. Including these factors, the general formula for the induced gluon number distribution can finally 
be written as 

x 



n / / m \ / m ^ 

-2 C (l, -,n) ■ ^2 B (m+l,-,n)(m,-,n) cos ^2^(k,-,n)^k -COS ^ ^(fe,-,w) & z h 



, k=2 / \fc=l 



(113) 

where Yla = is understood; 

We emphasize that Eq. (113) is not restricted to uncorrelated geometries as in [ ||. Also ^ allows 
the inclusion of finite kinematic boundaries on the q.; as well as different functional forms and elastic cross 
sections cr g (i) along the eikonal path. 
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For particular models of the target geometry one can proceed further analytically. For a sharp rectangular 
geometry, the average over the longitudinal target profile with 



• z n) — ~rz 0(L - z n )0{z n - Zn-i) ■ ■ ■ 0(z 2 - Z\) 



(114) 



leads to an oscillatory pattern [ [54| that is an artifact of the assumed sharp edges. 

A somewhat more realistic model [ |27], [2^] utilizes normalized exponential longitudinal distributions of 
scattering center separations: 



(115) 



This converts the oscillating formation physics factors in Eq. (113) into simple Lorentzian factors 

1 



/ 



dpcos ( ^a; (fej ... in )A2fe J =Rc Yl j 



^J- 1 + iw (fcr .. )n )L e (n) 



(116) 



In order to fix L e (n), we require that (z^ — z ) = kL/(n+ 1) for both geometries (114.115). This constrains 
L e (n)=i/(n + l) [11 1|. 

Note finally that the ratio of the medium induced to the medium independent factorization gluon distri- 
butions vanish at both small and large |k| 



k|=0 and oo dN^ 







(117) 



As |k| 0, dNW -► oo while all but one term in dN^ is finite. The potentially singular term cx 
(1/k 2 ) J Y^, j k • however vanishes due to azimuthal integrations. In the |k| — > oo limit, on the other 
hand, all Gunion-Bertsch currents B(..w..) — > 0(l/k 4 ) and hence vanis h fa ster than 1/k 2 from dN^°\ This 
observation partly justifies neglecting the kinematical |k| boundaries (120) in analytical approximations. 



V. NUMERICAL ESTIMATES 



A. Angular distributions 



Figs. 4a-c illustrate Eq. (113) for a exponential geometry (115) equivalent to a box of thickness L = 5 
frn with X g = 1 fm including kinematical constraints appropriate for a 50 GeV jet. We used adaptive 
Monte-Carlo integration [ |lj to integrate over the momentum transfers qi • • • q„. 

The medium induced gluon differential distri buti ons up to third order in opacity are plotted divided by the 
zeroth order in opacity, hard distribution Eq. ( |119[ ) for radiated gluons with light-cone momentum fractions 
x = 0.05 , 0.2, 0.5. Note that dN^ md ^/dN^ tends to vanish at both small and large k in accordance with 



Eq. (117). As a; increases the relative magnitude of the medium induced contribution decrease and the 
overall magnitude is set already by the first order in opacity result. Higher orders redistribute the moderate 
k 2 /fi 2 ~ 10 toward higher values due to elastic rescattering of the gluon, but they also tend to fill in the low 
k region with additional soft radiation. The overall angular pattern appears somewhat complex, but one 
must keep in mind that the variation is rather slow since the logarithmic scale varies over several orders of 
magnitude. 

Fig. 4d shows the actual rather modest size of the induced radiation contribution on top of the hard 1/k 2 
distribution for the case of x = 0.05. 



B. Intensity distribution and energy loss to first order in opacity. Analytical approach 
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1. Zero Opacity Limit 
The jet distribution in the absence of final state interactions is given by [ 

^^ (0) ®2^- d «' J ('P''P)' 2 ^W ' (118) 

In the Leading Pole Approximation (LPA) approximation [ [53) , the radiation distribution accompanying the 
such hard processes for a spin \ jet is given by 

dN^ _ C R a s ( x 2 \ 1 

X tedV™— \ l - X+ ^) V ' (119) 

where in the eikonal approximation the light-cone momentum fraction x = k + / E + « uj/E and k is assumed 
to be small compared to u>. For other spin jets, another suitable splitting function replaces the x dependence 
above. 
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FIG 4. (a)-(c) The medium-induced double differential gluon distributions are plotted vs. k 2 //i 2 up to 
first (diV (1) ), second (cW (1+2) ) and third (diV (1+2+3) ) order in the opacity (L/X g ) expansion, divided by 
the medium independent hard radiation (dN^) oc 1/k 2 . Curves are obtained numerically for exponential 
geometry with L/X g = 5, Ej et = 50 GeV and /j, — 0.5 GeV. The first three figures (a)-(c) are for typical 
soft, semi-soft and hard gluons respectively, x — 0.05, 0.2, 0.5. (d) The full gluon differential distribution 
up to third order in opacity dN m + dN {1) + diV (2) + djV (3) is shown for x = 0.05 for L = 0, 3, 6 fm. 



We consider radiation outside a cone with |k| > //. and with the upper |k| bound determined from the 
three body (jet+jet+glue) kinematics. The gluon kinematic boundaries are therefore 

k 2 min = /i 2 , k 2 ^ = nrfn[4^a; 2 ,4£ 2 a:(l-a:)] . (120) 

The radiation intensity integrated over this range of k is 



where |k| max and |k| m ; n are given by Eq. ( pL20[) . Note that the differential intensity is roughly uniform with 
the exception of the kinematical edges x — > and x ~ * 1 as seen in Fig. 5a. In the Leading Log of the 
Leading Pole Approximation (LLA) the radiative total energy loss of a quark jet originating from a hard 
vertex outside a cone defined by Eq. (120) is given by 

AE ( 0) = 4Cna lE E 
3n \x 

While this overestimates the radiative energy loss in the vacuum (self-quenching), it is important to note 
that AE^ jE ~ 50% is typically rather large. As shown below, the medium induced energy loss is small by 
compar ison . For a gluon jet one has to replace the quark splitting function q — > qg by the gluon one g — ■> gg 
in Eq. ( |m|) and use C A = N c « 2C F . 

The total energy loss should reduce to the medium independent one in the limit of vanishing opacity 

lim m = 1 (123) 
and for jets of asymptotically high energies due to factorization 

AE (tot) 

lim — — T^r = 1 ■ (124) 



2. First Order in Opacity Correction 



The first o rder (n = 1) in opacity contribution, dl^ /dx, to the induced radiation intensity can be read 
of from Eq. ( 113 ). The longitudinal coordinate average over the equival ent e xponential target profile (|l!5| ) 
is done with L e (2) — L/2. Including the q — > qg splitting function as in (119) we have 



dl^ C R a s 



dx 





x 2 \ 

















dk 2 
~k^ 



d 2 qi 



2k-q 1 (k-q 1 ) 2 L 2 
7r(q 2 + m 2 ) 2 \§x 2 E 2 + (k - qi) 4 L 2 



(125) 



To obtain a simple a nalyt ic result, we ignore the kinematic boundaries and set |k| m ; n = 0, |k| max = oo that 
is motivated by Eq. (117)). We also set |qi| ma x = 00 (i- e - M^// = M 2 )- This allows us to change variables 
q ' ee k - qi in Eq. fll25|) 



25 



diil} c ^(i- x+ x -\^e q,2L 



dx ir \ 2 J \ g J ^ l6x 2 E 2 + q' 4 L 2 

dk 2 f 27r d(j) 2k • (k + q') 



k 2 J 27r(q' 2 +k 2 + 2|q'||k|cos(0)+Ai 2 ) 2 
and express the integrand in the azimuthal (f> integral as a partial derivative with respect to k 2 
2n ^ 1 1 



(126) 



/ 2vr (q' 2 +k 2 + 2|q'| |k|cos(0) +A* 2 ) "^((k 2 + q' 2 + /z 2 ) 2 - 4k 2 q'< 

The remaining q ' integral 

d7« C R a s x 2 \ L _^ f°° , , 2 2 M 2 q ,2 L 



dx ir V 1 X+ 2)x g E J Q dq ' q' 2 +n 2 16x 2 E 2 + q<iL 2 ■ (12?) 

can be performed then analytically, resulting in 

dl& C R a s ( x 2 \ L 

'l- x+ \E — /( 7 ), 7=731?' ( 128 ) 



da; 7r \ 2 J X g AxE 



where 7 is a measure of the formation probability. The formation function /( 7 ) is given by 

7j^+27fog7) f 7T 7 if 7 «1 

A7J (1+7 2 ) \ 2 log 7 if 7 >1 • 1 yj 



It is the 7< 1 limit of the formation function Eq. (129) in which the the first order in opacity medium- 
induced intensity distribution reduces to a simple form with a characteristic quadratic dependence on L: 

dlW „C R a, l-x+£ L*y? 

— - - . l-L^UJ 

ax A x X g 

This formula breaks down at both x — > and x — > 1 because |k| max and |qi| max cannot be approximated by 
00 and because the small x approximations used above break down as x — ► 1. 

The induced radiative energy loss to first order in opacity in the framework of the above approximations 
is then given by 

A£(1) = C^LV log £ 

4 Xg fl 

This equation is directly comparable to Eq.(??) of Ref. [ [l2||. 



C. Induced radiation intensity to higher orders in opacity. Numerical results 

Whereas Eq. ( |131| ) displays analytically the main qualitative features of non-abelian energy loss, in practice 
at the finite energies available kinematical bounds do affect quantitatively the results. This naturally leads 
to a reduction relative to the analytic estimate. 

In Fig. 5a the induced intensity distribution dl^ md ^ /dx (sets of dashed curves) is compared to the medium 
independent dl^°\ We consider the example of a 50 GeV quark jet (fi = 0.5 GeV) at RHIC. The hard 
radiation intensity result is roughly constant with x, whereas for most of the x region the induced radiation 
intensity falls like l/x a , a ~ 1. We note that for relatively thin plasmas L/X g < 3 the medium induced 
energy loss remains small compared to dl^ . 
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FIG 5. (a) The radiation intensity distribution is plotted vs. the light-cone momentum fraction x of 
the gluon. We consider a 50 GeV quark jet in a plasma with screening scale fi = 0.5 GeV and X g — 1 
fm. The solid curve shows the dominant medium-independent radiation intensity. The medium-induced 
gluon spectrum is plotted for up to third order in opacity (dl^\ d/' 1+2 ' and d/^ 1+2+3 ' ) for opacities 
L/Xg = 1, 2, 3. (b) The absolute value of the orders in opacity dl^\ dl^ and dl^ that contribute in part 
(a) are plotted for the same energy and opacity L/X g —3. 



It is surprising how much the first order result dominates the induced intensity distribution [ Eq|. Fig. 5b 
illustrates how rapidly the contributions from higher orders in opacity decrease. While the differential 
angular distribution continue to show some weak sensitivity to higher order contributions in opacity the angle 
integrated intensity is much less affected beyond first order. The third order contribution is almost buried in 
the "numerical noise". For opacities L/\ g ~ 6 some more pronounced probability redistribution can be seen 
in the double differential level Fig. 4. However, when integrated over k 2 and then over x m \ n < x < x max to 
obtain the induced non-abelian energy loss AE (ind) , the first order result dominates very strongly as seen 
in Fig. 6. 
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FIG 6. The radiated energy loss of a quark jet with energy E jet = 5, 50, 500 GeV (at SPS, RHIC, LHC) is 
plotted as a function of the opacity L/X g . (X g = 1 fm, /j, — 0.5 GeV). Solid curves show the first order in 
opacity results. The dashed curves show results up to second order in opacity, and two third order results 
are shown by solid triangles for SPS energies. 



The induced energy loss shown in Fig. 6 is for quark jets with energies Ej et = 5, 50, 500 GeV typical for 
SPS, RHIC and LHC. Higher orders in opacity (L/X g ) n , n > 2 have little effect except at very low SPS 
energies. The kinematic bounds at SPS suppress very much the energy loss in comparison to RHIC and 
LHC energies. An analysis of the slopes as a function of the opacity, L/X g , shows that AE(ind) oc i 2±01 
at all energies even with finite kinematic bound arie s included. As a measure of the deviation of the simple 



first order analytic estimate, we generalize Eq. (131) as follows 

-2„2 



N(E) Xg fl 



If the kinematic bounds are ignored as in Eq. ( J130| ), then N(E) — 4. Including finite kinematic con- 
straints cause N(E) to deviate considerably from this asymptotic value. We find numerically that 
N(E) = 7.3, 10.1, 24.4 for E = 500, 50, 5 GeV. Together with the logarithmic dependence of energy, these 
kinematic effects suppress greatly the energy loss at lower (SPS) energies as seen in Fig. 6. This is in contrast 
to the asymptotic, energy independent result quoted in Ref. [ Q| where the finite kinematic bounds were 
neglected. Our numerical results, however, agree with Ref. [ Ofnear LHC energies as reported in [Eg. 
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APPENDIX A: DIAGRAMS Mi, ,o = G Xi, , Mi,i, = GoG^Xi^Gi AND Mi, ,i = G Xi,i 

To make contact with the results in Ref. [ |27], we present explicit calculation of the first nontrivial 
diagrams shown in Fig. 7. 

M lj0;0 e k,c M 1>1>0 k,c M lf0 ,i 

Sqi' a i Sqi' a i fqi.ai ' c 



FIG 7. Three "direct" terms Mi,o,o = GoXi j0 , M M , = GoG _1 Xi,oGi and Afi,o,i = G Xi,i contribute 
to the soft gluon radiation amplitude to first order in opacity L/X oc <j e i/A±. 



As a first application, consider the one rescattering amplitude Mi,o,i, in the notation of Ref. [ 27 



M li0 ,i= / ^ I iJ(p + k~q 1 )e^ +k -^A 1 (p,k,q 1 )V(q 1 )e^ X 
(2tt) 4 



x iA(p + k- qi)(-i)A(k - q x ) 



« J(p + fc)e l (f +fc ^° [c, a x ]T 0l (-») j ^ e -*n-(*i-*<0 2 g a e • (k - qi ) x 

x E+ J ^v(q lz ,< ll )A(p+k-q 1 )A(k-q 1 )e~^^-^ ■ (Al) 

The longitudinal momentum transfer integral 

h(p,k^ uZl -z )= J ^v( qiz ,c il )A(p + k-q 1 )A(k-q 1 )e' 1 ^^-^ (A2) 

can be performed via closing the contour below the real axis since Z\ > z . In addition to the potential 
singularities at ±«^i (/if = = qf + M 2 ), the two propagators have four poles in the q\ z plane, which up 
to leading correction are located at 



Qi = E + + ie , q 2 = -u> - ie , 

q~3 = k + — uji + ie , q 4 = -cj + lji — ie . 



(A3) 



The residues give 



while 



e iu! (zi-zo) e i(uJo-i*Ji)(zi—zo) 

Res(g 2 ) ~ -v(-u ,<il) E+k+UJi ' Res (94) ~ v(wi - w 0) qi) E +j t + U j 1 ' ( A4 ) 

Ana e~ wO* 1- z °) 

Res(-ifix) w ■ \2 ' ( A5 ) 

(— 2i[ii)£j + k + (— z fix) 

where we assumed that k + ^> /ii ^> wj. In the well-separated case where \i{z\ — Zq) = /iA ^> 1, this potential 
residue is exponentially suppressed and therefore 

e iuj a (zi-z Q ) . , 

h(p,k,qi,zi - zq) ps i_ ^-j— ^(-w ,qi) - - w ,qi)e"^ l(zi " 2n) J 

p iu (zi-z ) . N 

-^^w^-^-^ 1 ^ 1 "^) • (A6) 
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We thus recover in this eikonal limit the time ordered perturbation result in [ ^7], £6[ 

M w = J(p)e ipX ° H) / |^(0,qi)e- iqi - bl x 

x 2ig s - ' (k ~ ^ e<(«o-<-0*i( e *-i«i _ e *-i*o) [ C)ai ]r ttl , (A7) 
(k qi ) 

where bi = xi — xo. Similarly, we recover the other two amplitudes for one center, e.g., 

M ljlfi = J[p)e ipXo (-i) J ^L v (0, qi )e~^ x 

x (-2* 5s ) ^ e 1 ^ Cai T ai . (A8) 

Note that while the direct (Mi )m] ;M| m , ;/ ) ensemble averages lead to non-vanishing modifications to the 
initial hard spectrum, |Mo| 2 cx |J(p)| 2 /k 2 , the interference between M$ and Mi,/ i?T1 vanishes because the 
ensemble average is over a color neutral medium with TrT Q (i) = 0. 

APPENDIX B: DIAGRAM M 2 ,o,3 = G X hl X 2 ,i 

Consider next the gluon two-scatterings amplitude M 2 $^ in [ [26fl . Fig. 8 shows that for inclusive processes 
two interesting cases arise. 



^2,0,3 M 2; o,3 

_^ZjQ Z Z^ Z2 ^xjZq Z Z^ Zj 

•-v - •-^ + 

q b ai® ® q 2 ,a 2 qi> a i ^ lilB ^ A2&2 

FIG 8. A/2,0,3 = GoX lt 1X2,1 "direct" contributes to second order in opacity cx (a e i/A±)' 2 , whereas M£, ,3 
A/2, 0,3(^2 = 21) = GoOi,i "contact-limit" may contribute to first order in opacity oc (a e i/A±) 1 as well. 

In the Feynman diagram approach 

7^77^7 iJ(p + k - gi - q 2 )e^+ k -^-^V{ qi y^V(q 2 y^ x 

x A 12 (p,k,qi,q 2 ) iA(p + k - q 1 - q 2 )(-i)A(k - q x - q 2 )(-i)A(k - q 2 ) x 
^J(p + k)e^ +k ^ [[c,a 3 ],ai](T 0a (2)T 0l (l)) x 

x H) / |^H) /^2 t5se .(k- qi - q2 )e-- b le --^ x 

d<?i* dq 2z (E + k+)v(q 1 )v{q 2 )e- i i 1 ' (zi - Z0 h- lq2 ^ Z2 - z ^ 



(2tt) (2tt) ((p + fc - ft - q 2 ) 2 + ie){{k q 2 f + ie){{k - q 2 f + ie) 



(BI) 



where = Xj — xo arc transverse impact parameters, and we used the soft glue and rescattering kinematical 
simplifications Eqs. ( p9|j3^ ), e.g. J(p + k — q\ — q 2 ) w J(p + k) w J(p). For the qi 2 integral, it is convenient 
to rewrite the phase as 
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The first longitudinal integral is closely related to Eq. (JA2J) 



I 2 (p, fc,qi,q 2 ,2i - zq) 



2tt ((p + k — qi — q 2 ) 2 + ie){{k -q x - q 2 ) 2 + ie) 



(B2) 



Since z\ — zq we again close the contour in the lower half q\ z plane a nd n eglect the pole at —ifii. 

The remaining q\ z poles are shifted by —q 2z and qi — > qi + q 2 relative to Eq. (A3): 



qi = -q 2z + E + +ie , q 2 
q~3 = -q 2 z + fc+ - w (12 ) + ie , q 4 



-qiz -wo-te, 

-qiz - +W(i2) - *e , 



where now fc + U( 12 ) = (k — qi — q 2 ) 2 . The residues at q 2 , q^ then give 

p iuio{zi — Zo) 



(v(-q 2z - cj , qi) - »(-&, - + W(i2), qi)e~ <W(I2)(zi_Zo) ) 



(B3) 



(B4) 



where we have neglected C(exp(— /iA)) contributions. This differs from Eq. ( A6) mainly in that the potential 
is evaluated near —q 2z , which still remains to be integrated over, and lo\ — > wm 2 )- 
Next we need the following critical q 2z integral 



7 3 (fc,qi,q 2 ,z 2 - z\) 



dq 2z v(-q 2z + 5uj, qi)v(q 2z ,q 2 )e 



(B5) 



2tt ((fc-g 2 ) 2 +«e) 
Fortunately, we are interested in only two extreme limits: 

• The limit of well-separated scattering centers z 2 — z\ 3> 1/fi ; 

• The special "contact" z 2 = z\ limit to compute unitary contributions. 

In the first case, we can proceed ignoring the q 2z = —ifix and — ifi 2 potential singularities. Then only the 
pole at q 2z = lu 2 — loq — ie contributes and yields 



I 3 (k, qi , q 2 , z 2 - z x » 1/n) » -pr«(0, qi)«(0, q2 ) e -^—°)(^i) 



(B6) 



In this case we recover the result quoted in Ref. [ |2^] for M 2i o,3 

d 2 qi 



M 2 , ,3« J(p)e""°(-i) 



z;(0, qi )e-^- bl H) 



rf 2 qi 
(2tt) 2 



•y(0,q 2 )e 



— iq2 b 2 



x (2i&) 



e- (k-qi -q 2 ) 
(k- q x - q 2 ) 2 

x [[c,a 2 ],a 1 ](T a2 T ai ) . 



(B7) 



In the general case (including the special contact case with z 2 = Z\) both q 2z — — i/i 2 , —ifJ-i singularities 
in the Yukawa potential contribute together with the pole at q 2z — lu 2 — luq — ie, resulting in 



7 3 (fc,qi,q 2 ,z 2 - zi) ~ ( v(0, qi)«(0, q 2 ) 



-i(0)2— U)o)(Z2 — zx) 



'2(m?-m!) 



3-^2(^2-^1) 



j-/il(z2-Zl)g-i<5w(^ 2 --2 : l) 

/2 



(B8) 



For z 2 — zi = A>l//i this reduces to Eq. ( B6 ) . For the special contact contribution z 2 — z\ = it reduces 
to 



J 3 (fc,qi,q 2 ,0) 



2fc+ 



w(0, qi)u(0, q 2 ) , 



(B9) 



i.e., exactly \ of the strength in Eq. (B6). The contact limit of this amplitude is therefore 



31 



M c 2 ^^J( P y^{-i) J ^ u( 0, qi )e-^' bl H) J ^«(0,q 2 )e- iq2 - b2 



x 



1 ( 2i ^) ^ V elUJ ° Zl (l-e^ (12)(21 " Zo) ) [[c,a 2 ], ai ](T Q2 T ai ) . (BIO) 



(k - q : - q 2 ) 



After this explicit derivation of the factor i in the contact limit, we can generalize it to any functional 
form of the potential as follows. First we assume that a contact interactions only contributes for q*x = 
— q 2 , corresponding to no net transverse momentum exchange inside the potential function only, i.e. 
v(qi z ,qi) = v(— q 2zi — q 2 )- However the distinction has to be kept in the pro paga tors and phases. In this 
case the potential function f(qi), does not appear in the first integral I 2 , Eq. (B2), which is modified to 



h{p, fc,qi,q 2 ,zi - zq) = J 



e 

1 



dq lz e -i(«ii,+<»,)Oi-*o) 
~2^T ((p + k - qi - q 2 ) 2 + ie)((k - ~q\ - q 2 ) 2 + ie) 

iuio(zi-zo) 



1 _ e -<"(ia)(*i-*o) \ _ ( B11 ) 



S+fc+CJ(i2) 

This is independent of the functional form of the potential. The second integral, I3, then is modified to 

fdq 2z \v{q 2zi q 2 )\ 2 e'^^-^ 



h(k,qi = q 2 ,z2 - ^1) = / 



2ir {{k-q 2 ) 2 +ie) 
1 fdq 2z \v(q 2zi q 2 )\ 2 e-^^-^ 



k + J 2tt q 2z — Suj + ie 



(B12) 



where z 2 — z\ = A and we have used the finite range of \v(q z , q)\ 2 , i.e. v is limited within (— /i, to 
neglect the residue at q 2z ~ k + . For large z 2 — z\ compared to the range of v, we close below and obtain the 
dominant residue 



h~-^H0,q 2 )\ 2 e-^-^ s " (B13) 

independent approximately of the actual form of v as long as all singularities of v in the lower half plane 
have imaginary parts — ifa with /ijA = [ii{z 2 — Z\) 3> 1. 

To extract the contact limit, z 2 — z\ = 0, we can integrate along the real q 2z axis and use the Dirac 
representation of a pole approaching the real axis: 

1 q 2z -Suj 

nrd(q 2z - buj) .. 



q 2z - Suj + ie (q 2z - Suj) 2 + e 2 
This gives 

^-)l3,-l Kfciq2 )P + ^ |i|, fe + ^,q 2 )P ^ C ° S fe A )-^ SinfeA) , (B14) 
2 J^^ 2ir q 2 z + e 2 

where we have introduced q z — q 2z — Suj. For our high energy eikonal applications Suj <§; fx, and therefore we 
can now expand \v\ 2 around Suj = 0. The correction of 0(6u)/n) can be neglected. This leads to a vanishing 
real part since that part of the integrand is odd. On the other hand, the imaginary part has a finite e = 
limit 



d( lz, 2 sm(g z A /O if /iA = fx(z 2 — Zi) —> Q 

-^-Hqz,q2)\ = = i i , , n \, 2 -t ^ , ^ > B15 ) 

, 2tt q z { -l\v(0,q 2 )\ 2 if ^A = fi(z 2 - Zi) -> 00 

where fi is the range of the potential v as indicated above. We see that while the detailed interpolation 
between the asymptotic limits depends on the actual functional form of v, the limits Eqs. (B6 B9) and the 
factor of i reduction in the contact (A = 0) limit is very general. 

We thus confirm that O^i diagrams where the two legs are attached to the gluon line can be obtained 
from the higher order direct diagram by setting z%+i = Zi and reducing the strength twice (i.e. times |). 
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APPENDIX C: DIAGRAMS M 2 , ,o = G a X lfi X 2 ,a AND M 2 , 2 ,o = G Xi,oG" 1 X 2 ,oG2 



In those graphs it is the jet rather than the gluon that suffers two sequential scatterings as seen from 
Fig. 9. 



k,C 



M 2,0,0 ^ 
W!q A^Zj Z 2 

— 1- 



z I 



qi.ai 



q2^ a 2 



® — ^-r r- 



qi»ai 




q2,a 2 



k,c 

M 2 , 2 ,o ^ M^o 
® 1 8—^ 



s I z P 




qi,a!® § q 2 ,a 2 



FIG 9. M2,o,o = Go-Xi, 0X2,0 and M2,2,o = Go-Xi,oG 1 X2,oG2 graphs in the well-separated case together 
with their z 2 = zi limits M^o = GoOi.o, A/| i2 ,o = GoG _1 Oi,oGi. 

It is straightforward to write the expression for the amplitude 

M 2 , ,o - J ^4^4 *J(P + k - 91 - q2)e^ +k -^-^ X0 V( qi )e^V(q2)e 1 ^ x 

x (-iE+) 2 ig s (2p + k)^ iA(p + k-q 1 -q 2 ) iA(p - qi - q 2 ) iA(p - q 2 ) x 



x a 2 a lC (T a2 T ai ) 

r) 2 



^)^°H)/0 e - iqi - bl H)/^e — < 



2i gtl (e-k) eiuolQ ^^^J {E+) 2 x 

a; 

dfe dq 2z w(qi)w(q 2 )e-^( Zl - Z0 )e-^( Z2 - Z( ') 



J (27r)(27r) 



r) + fc - gi - q 2 ) 2 + «e)((p- <7i - <72) 2 + «e)((p - <72) 2 + ie) 
In this case we define 



Since z\ — z ^> 1/^, we neglect the pole at — iyUi. The remaining q lz poles are 

9i = ~1iz + E + + k+ + ie , q 2 = -<? 2z - o> - ie , 
93 = -fe + E+ + ie , <?4 = -q 2z - «e , 



(CI) 



, dg lz < glz , qi )e-^+^>^- z "> 

J 2 (p,A:,qi,q 2 ,zi - z )= / -— — 77 , . 777 x 2 , . ; •• (C2) 

2?r {(p + k - q x - q 2 ) 2 + ie)((p - q x - q 2 ) 2 + ie) 



(C3) 
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where we discarded (p + k — qi — q 2 ) 2 /E + relative to ujq. The q 2 , qn residues then give 



.«(-«!«, qi) 



o iui (z 1 -z <l ) 



■ 



(C4) 



(£+) 2 ^> V 

Not that u>o has been neglected in the potential relative to fii. The second integral, J3, is then modified to 

dq 2z v(-q 2z ,cn) v(q 2z ,( l2 ) e -W*a-*0 



^(p, qi,q2,22 - zi) 



(27T) 



u(0,qi)w(0,q 2 ) x 



((p-q 2 ) 2 +ie) 

1 if (LtA = /i(z2 — ^1) — > 00 



if /iA = /i(z 2 — z\) 







(C5) 



Note that E + ojq = k 2 /x. With the help of Eqs. (C4,C5) in the case of well-separated scattering centers we 
obtain 



M 2A o = J(p)e^°H) | |^(0, qi )e- iqi bl H) | |^(0,q 2 ) 



(27T) 

2», gs (e-k) 4 



e ,q2 b2 x 



') a2aic(T a2 T 0l ) 



(C6) 



We thus recover the result from [ p7| . In the contact limit there is extra factor of | relative to the naive 
expectation and Eq. (C6). Analogous calculation leads to the expected result for M 2 2 . 

We thus arrive at a conclusion similar to the one in Appendix B, i.e. in the case when the two legs of the 
virtual contribution are attached to the quark line, O^O: there is an additional factor of | relative to just 
setting z i+ i — Zi in the higher order well-separated case amplitudes. 



APPENDIX D: DIAGRAMS M 2 , ,i = G X ltl X 2 ,o AND A/2,0,2 = G X lt0 X 2 ,i 

In the case when one of the hits is on the parent parton and the other hit is on the radiated gluon we 
encounter a different situation. Explicit calculation shows this on the example of Af| 1 = M 2] o,i(^ 2 = Zi) 
in Fig. 10. 



M 



2,0,1 

Z Z Zj 



Zl 

-§- 



i p 

% I k,c 



qi.ai' 



q2,a 2 



M C 2 ,o,2 



J-0 z z l 



l-rrmmmm-n l. - 



FIG 10. M^o.i = GoOi,2 and Af£o,2 = GoOi,2 topologically indistinct contact diagrams. There are no 
additional factors of = arising from the integration in taking the Z2 = Zi limit in A/2, 0,1 = G0X1, 1X2,0 an( l 
Af 2 ,o,2 = G0X1, 0X2,1 . 

M 2 , .i = / ^^TT + * - 91 - g2 )e l(p+ ^ 9l - 92) ^^( 9 i) e ' i9ia:i F(g 2 )e^ a;2 x 
x (-iE+)A 1 iA(p + k-qi- q 2 ) (-i)A(fc - qi)iA(p - q 2 ) x 



J(p)e^°(-i) j ^e-^(-i) J 2 x 

x 2ig s (e- (k-qi)) e luJ ° z ° a 2 [c, ai ]{T a2 T ai ) (E+f x 



dqiz dq 2z v^v^e'^ 1 ^ 1 " 2 ^ e~ iq2z(z2 ~ za) 

(2tt) (2tt) (( P + k- qi - q 2 f + ie){{k - qi f + ie)((p - q 2 ) 2 + ie) 

We perform the q\ z integral first 



(Dl) 



h p,fc,qi,q 2 ,zi - z ) = / -5— ,, , , x 2 . ■ ,,,, WTT\ ' D2 > 

7 2?r + k - qi - q 2 y + ie)((k - qi) z + ie) 

The pole at —i[i\ is again exponentially suppressed. The poles of interest in the lower half plane are 
qi z = —q 2z — ujq — ie and q\ z — —ujq + lo\ — ie. Taking the residues leaves us with 

iwo(zi-zo) 

72 = 1 — ~M-^ - w o, Qi) - - Wo , qi )e-^+^ Zl -^) . (D3) 

E+k+(q 2z +uji) 

It is important to notice that there is no pole at q 2z — —uji in Eq. ( |D3| ). The remaining integral over q 2z is 

f dq 2z 1 / er iq2 ^^ Z2 ~ Zl ^ 

h{p,k,qx,c^,zx - z 0) z 2 - z x ) = / — — ■ . »(-g 2a ~ ^o, <l 1 )v(q 2z , q. 2 ) 

J 2tt q 2z +LUi \(p-q 2 y+ie 

e -i(q2z(z2~Zo)+Ul(z 1 ~ Zo )) \ 

7 Try v(ui-uJo,q 1 )v{q 2z ,q 2 )) . (D4) 

Unlike in previous examples, we now show that no factor of \ arises in the contact limit z 2 = z\ > Zq. The 
poles in the lower half plane are q 2z = — ie, q 2z — —i(i2, and q 2z = — ojq — ifi 2 and the corresponding residues 
contribute 



Res( _ ^ r,(0 1 q 1 M0,q a ) (1 _ ^ iu>1 ^ Zo)) 



Res(— ifi, 2 ) 



E+nl(2i(i 2 )(nl - nl - 2iuj fj, 2 ) 



U vas f e -^(z 2 - Zl ) e +i^( Z2 - Zl ) 

The exponentially suppressed contributions oc exp[— /i2(^i — Zo)] in the second and third residues have been 
neglected, and we made use of the approximation u>i -C fij. In the limit of well separated scattering centers, 
the second and third residues are also exponentially suppressed. In the contact limit, on the other hand, both 
the second and third residues in Eq. ( p5| ) contribute. However, impact parameter averaging via Eq. ( JTgj ) 
sets \i\ = fj, 2 , and the those two contributions cancel exactly. The result is then 

M| A1 =J(p)e^H) J |^L e -^ b M0,qi)H) J |^|e-^ b M0,q 2 ) x 

x 2ig s 6 ' (k ~ ^ e <("o-"i)»i(e<"i'i - e ^ z «) a 2 [c, a{\{T a2 T ai ) . (D6) 

We emphasize that in this case, there is no factor of ^ emerging in the z 2 = Z\ contact limit Oi 2 when one 
of the legs is attached to the quark line and the other on to the gluon line. 

Similarly, the diagram M 2 .o,2 leads to the result quoted in Ref. [ ^?J. In the contact limit, it reduces to 
-^2,0.2(^2 = Zi) — M 2 fl.\(z 2 =21). However, it is important to point out that in the contact limit only one 
of the diagrams must be taken into account in order to avoid over counting. This can be directly seen from 
the expansion in the powers of the interaction Lagrangian. Alternatively, we get the correct contact answer 
by taking the limit of both cross contact diagrams and multiplying each by | ■ 



35 



APPENDIX E: ZERO MEASURE CONTACT LIMIT OF 

M a ,i,o = Xi, GiX 2 , AND M 2 ,i,i = Xi, GiX 2 ,i 



In calculating the different contributions coming from two interactions with the same potential centered 
around Xi we have to take into account the two graphs given in Fig. 11, where one of the hits occurs before 
the gluon emission vertex and the other one after. 




M' 



2,1,1 
z o 



qi,ai 




^ p 

4? _ k,c 



FIG 11. Diagrams M 2 c il = X lfi GiX 2 ,o(z2 = Zi) and M^ 1:1 Xi,oGiX 2 ,i(z2 = zi) of O(0) according to the 
time-ordered perturbation theory. 

In the framework of time-ordered perturbation theory from [ ^6) the graphs are identically zero in the 
contact limit because 1 dt ■ ■ ■ = 0. We here present a more detailed study of the validity of this argument. 
In the contact limit ^2,1,0 — > M 21 = X\ qG\Xi q. We will show that this contribution vanishes in the 
small x = k + /E + limit 

M 2 c jlj0 = J iJ(p + k - qi - q 2 )e^+ k -^-^V{ qi y^V{q 2 y^ x 

x (-iE+) 2 ig s (2(p - q 2 ) + jfe^e" iA(p + k- qi -q 2 ) *A(p + k - q 2 ) iA(p - q 2 ) x 
x a 2 cai(T a2 T ai ) 



2tt) 2 v 'J (2tt) 

x 2i 5s ^ e*** a 2Cai (T a2 T ai ) (£+) 2 | ^y^y K^M^) x 

e -i(gi^+<32z)(zi-zo) / 1 1 \ 

(p + k - qi - q 2 ) 2 + ie \ (p - q 2 ) 2 + ie (p + k - q 2 ) 2 + ie ) ^ 

In writing this we have made use the simplifying soft gluon and weak interaction approximations in Eq. ( p9[ ) . 

We can now perform the q\ z integral, which brings a factor of 1/E + and sets the overall phase to 0. We 
have to now perform 

h= [ 7T w(-92 2 ,qiM<?2 2 ,q 2 ) (7 ; —r 1 „ - . J • (E2) 

J 2?r \(p-q 2 ) 2 + ie [p + k — q 2 ) z + ie/ 

We close in the upper half plane where there are four poles q 2z — ifi-y, q 2z — ifi 2 , q 2z = E + + ie, q 2z — 
E + + k + +ie. The residues at the poles involving E + are suppressed by a factor (1/E + ) A and have negligible 
contribution in the high energy approximation. The remaining two residues give 

(47ra s ) 2 / k+ 



, . (4-7ra s ) 2 / k+ , 
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We have discarded uj relative to fa and carried the k + /E + expansion to first order. It is important to notice 
the sum of the residues is not singular for the special case of interest when the contact limit contribution is 
averaged over the transverse position of the scatterer, i.e. fa = fa>. The above contribution is suppressed by 
an 0(k + 1 E + ) factor relative to the other graphs. In the high energy E + — ► oo limit, where the time-ordered 
perturbation theory works, ^2,1,0(22 = z i) ~ 0. Similar calculation for M^^i — > M$ % xXx^GxXi^ shows 
suppression of the order 0(/j,/k + ) and is approximately zero according to the k + ^> /1 regime ( pq ) in which 
the problem is set up. 

This leads us to the conclusion that the naive "zero measure" argument for such graphs is approximately 
valid for our kinematics. 



APPENDIX F: AMPLITUDE ITERATION TECHNIQUE TO SECOND ORDER 



To illustrate the general iteration procedure and check the general iteration results using the reaction 
operator approach, we construct the classes for rank 1 and 2 explicitly. 

For the two rank 1 classes, we apply the operators D\ and V\ on the hard vertex amplitude Eq. ( |39| ) once 
to obtain 



GqDi = H e waza 01c - Ci e^o^i^w) ^ j _ Bl e **a [ C; 0l ] f 
G Vi = Cr + Ca h e luJ ° z ° C-—C-L e <(<"o*i-«i*io) c _ £A. Bj e iuJ ° Zl c 



(Fl) 
(F2) 



We recall that Zij = Zi — Zj, Azi = Zi — These are all amplitudes needed for the first order in opacity 

L/X calculation. 

Some of the rank two classes are obtained from the rank 1 classes Eq. (F2) through relabeling, i.e. 
G D 2 = G Di(l -> 2), G V 2 = G Vi(l -> 2). The rest are readily derived from Eqs. through our 

iteration scheme Eqs. ([76 J 



G0-D1-D2 = -He" 020 a 2 aic- Cie l(a,oZl - WlZlo) a 2 [c,ai] - C 2 e 1 ^" 22 '^ 2 ^ ai [c, a 2 ] 
-Bi e lW0Zl 02 [c, 01] - B 2 e iu » Z! [c, a 2 ]a x 



G D 1 V 2 



-C ( i 2) e^2 — a * 3 i-« (i 2 )^o) [ [Cj a2 ] ; ai ] _ b 2(12) e «<-o*2-<**») [[Cj a2 ] ; 0i] 



(F3) 



Cr + Ca 



Cie l(a)0Z1 ~" lzl0) [c, ai] + C 2 e 



i(wQZ2-W2Z2o) 



a 2 ai c, a 2 



+C (12) e <(<-o*.-**i-«Wio) a2 [[ C) a2]) ai] + B2(i2) e ih^^^) a2 [[ C; 0a]j 0l ] 



(F4) 



Bi e^ 1 a 2 c + ^ B 2 e 1 ^ 2 [c, a 2 ] 



C A 



C A 



-ir C (12) e 1 ^ 022 "" ***->W) [ c , fl2 ] - B 2(12) e l < 



(W0Z2- W222l) 



c, a 2 



(C« + c A f 



He 



(F5) 
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Ca{Ca + C R ) e i( UoZl -uiizi ) c + c a(Cr + C A ) c ^(^022-^2220) c 
4 4 



Ca(Cr + Ca) 



Bie 11 



c + 



CrCa 



B 2 e 11 



/^2 /-<2 

-_d c (12 ) e^" 022- " 2221 ""* 12 ) 210 ' c B 2 (i 



(12) e l(a,0Z2 -" 222l) c 



(F6) 



With this explicit construction of the relevant classes, we can compute the differential gluon probability up 
to second order in the opacity expansion. Similar calculations can be performed for the Gunion-Bcrtsch case 
as well. 

To first order in opacity L/X the probability to radiate a gluon from either quark or gluon jets is proportional 

to 



1 



Tr 



G D 1 D\Gl + 2Re(G V 1 Gl) ) = C R C A { [ (-2 C x • Bi) (1 - cos^A^)) ] >„ 

J / V 

= CrC a J dq x (« 2 ( qi ) - S 2 ( qi )) [ (-2 d ■ BO (1 - cosKAzO) ] , (F7) 

where the 17 terms in Eqs. (fj^Jiil ) collapse to a single term with color trivial factor CrCa<1r interpretable 
as gluon final state interaction only. We note that if the normalized potential v 2 (qi) is replaced by £> 2 (qi) — 
<5 2 (qi) the result of Eq. (F7) remains unchanged. 
To second order in opacity we find that 

i ( Tr [G DiD 2 DlDl& + (GoV^DlGl + h.c.) + {G q D x V 2 d\gI + h.c.) 

+(G V 1 V 2 Gl + h.c.) + {GaVtV^Gl + h.c. 
= C R C\ ( [2Ci ■ Bi (1 - cos(wiAzi)) + 2C 2 ■ B 2 (cos(w 2 Az 2 ) - cos(w 2 (Azi + Az 2 )) 

- 2 C(12) • B 2 (12) (l - COS(W(i2)A2i)) 

- 2 C( 12 ) ■ B 2 (cos(w 2 Az 2 ) - cos(w (12 ) Az 1 + w 2 Az 2 )) ] ) v 
= CrC\ H f d qi (v 2 ( qi ) - 6 2 ( qi )) [- 2 C (12) ■ B 2(12) (1 - cos(o; (12) A^)) 



i=l,2 1 



2 C ( i2) • B 2 (cos(w 2 Az 2 ) - cos(w (12 ) Az! + w 2 Az 2 )) ] 



(F8) 



Eqs. (F7 F8) have simple physical interpretation. We showed that for both quark and gluon jets a simple 
color trace cx C\ survives. In general this can be seen from Eqs. 



,|82j,|81|) . Thus the distribution of the 
radiated gluons is interpretable as the interference of the Cascade amplitude that has the knowledge of all 
final state interactions with the effective color currents, i.e. the Gunion-Bertsch amplitudes, originating 
at the scattering centers and also having all possible subsequent interactions with the gluon line. The 
formation physics effects arise from phase differences that store the information on the cumulative formation 
lengths before the effective color current emission. The averaging over the momentum transfers is seen to 
be effectively over a modified potential that vanishes on integration. 
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